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ABSTRACT  ■ V' 

The  advances  made  in  the  development  of  high  .speed  electronic  com- 
puting and  equipment  have  increased  the  importance  of  analytical  methods 
in  precision  photogramme  try. 

The  analytical  treatment  of  orientation  and  calibration  6f  a p’hotyj- 
grammetric  camera  is  discussed  from  the  point  of  view  o.f  the  geometrical 
configuration  as  well  as  the  least  squares  adjustment.  The  solution  can 
be  applied  to  measuring  methods  in  both  ground  and  aerial  photogrammetry, 
The  projective  relation  between  two  plur.fts  is  interpreted  as  a special 
case  of  the  general  problem,  which  defers  with  the  relation  between  the 
plane  of  the  plate  and  the  spatially  arranged  objects. 

The  solution  is  suitable  for  electronic  computing  devices  and  is 
marked  by  its  computing  economy. 
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TMTRiTU  RATION 


The  stereoscopic  evaluation  of  pairs  of  photographs  using  universal 
precision  plotting  instruments  had  become  an  established  technique.  These 
modern  photogramme  trie  methods  have  by  far  their  widest  application  in 
geodetic  work,  namely,  in  the  triangulation  of  control  points  and  in  the 
compilation  of  topographic  maps.  Besides,  in  many  scientific  institutes 
and  research  laboratories,  photogramme trie  methods  are  in  use  which  would 
provide  a good  approach  to  many  different  measuring  problems.  The  re- 
markable advance  in  the  development  of  powerful  precision  surveying  lenses 
with  relatively  wide  angle  of  view  and  the  introduction  of  intefciiangeable 
glass  pl.ates  of  highest  quality  of  flatness  a'nd  emulsion  have  made  the 
rjqdern  photogrammetric  camera  with  its  extreme  inner  stability  another 
tool  for  precision  measurements , The  increasing  use  of  large  scale  photo- 
graphs by  airborne  cameras  for  precision  surveying  and  similar  application 
of  such  cameras  on  the  ground  and  in  the  air  for  precision  measurements  in 
non-topographic  fields  have  made  it  necessary  to  study  the  problem  of  ori- 
entation of  a phot bgrammetrie  camera  and  the  propagation  of  residual 
errors  of  the  process  of  orientation*  The  error  theory  of  stereo-photo- 
gramme try  has  already  been  discussed  by  several  authors  [p  .g . (£),  (13)'~| 
and  valuable  detailed  information  on  the  subject  has  been  published 
j_e,g.  (1),.  (2),  (3)].*  However,  these  studies  are  limited  almost  ex- 
clusively to  the  problem  of  determining  the  exterior  orientation  of  two 
approximately  vertical  photographs  by  means  of  a universal  plotting 
instrument.'  This  specialization  is  consistent  with  the  importance  of 
this  problem  in  the  field  of  topographic  photogrammetry „ However,  there'1 
is  an  indication  that  in  the  future  theoretical  photogrammetry  must  treat 
its  problem  not  only  in  connection  with  certain  types  of  evaluation  ma- 
chinery, but  must  continue  analytical  treatment  where  in  the  late  twenties 
the  optical  mechanical  reduction  methods  took  over.  At  that  time  the 
numerical  evaluation  seemed  too  laborious  to  be  considered  for  practical 
purposes . Today  the  development  of  an  economical  numerical  solution  of 
certain  photogrammetric  problems  appears  to  be  within  reach  with  electronic 
computing.de  vices  • 

The  continuously  increasing  precision  in  the  elements  of  photo- 
graimuetric  cameras  has  made  the  camera  the  most  precise  component  in  the 
chain  of  steps  leading  from  the  photograph  to  the  coordinate  determina- 
tion. Consequently,  evaluation  instrumentation  must  enable  ua  to  measure 
the  plate  coordinates  of  certain  points  - either  for  the  purpose  of  ori-  « 
entation  or  final  coordinate  determination  - with  corresponding  high 
accuracy.  It  is  well  known  that  the  monocular  or  stereoscopic  comparator 
designed  according  tc  Abbe's  comparator  principle,  and  used  with  pre- 
cision grid  plates  for  calibration,  is  an  adequate  plate  reading  device  to 
meet  precision \reqirLrement3 , However,  the  use  of  such  an  instrument  balls 
for  computation  of  the  desired  results  from  the  measured  plate  coordinates 
by  rather  complex,  formulas  and  at  high  speed.  In  addition  to  the  re- 
quired automatic  high  speed  computer  an  automatic  comparator  reading  de- 
vice, leading  to  an  automatic  computer  input  and  an  automatic  data  out- 
put, is  required.  In  view  of  the  possible  development  of  small  electronic 
computers  together  with  other  electronic  recording  deyices,  these  require- 
ments may  possibly  be  met  at  a cost  comparable  to  that  of  present  day 


•*»  Bibliography  at  end  of  paper. 
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universal  plotting  -machines*  It  is  with  these  thoughts  in  mind  that  the 
aimly  fclcal  determination  of  the  orientation  of  a photogrammetric  camera  i,3 
treated  in  detail,  since  this  step  is  the  common  basis  for  all  photo- 
gramraetric  measuring  methods.  In  addition  to  the  use  of  independent  photo- 
graphs for' special  problems,'  the  stereoscopic  method  may  be  considered  as 
a combination  of  two  independent  cameras.  The  geometrical  conditions  as 
they  exist  between  tine  two  bundles  of  projective  rays  can  be  introduced 
analytically,  by  adding  such  condition  equations  in  the  process  of  the 
analytical  reduction  as  are  given  by  the  functions  expressing  the  relative 
orientation.  Such  an  approach  seems  to  be  useful  for  analytical  treat- 
ment and  shows  again  the  importance  of  the  basic  problem,  namely,  the 
Orientation  of  the  individual  photogrammetric  camera. 


it 


THE  PHOTOGRAMME TRIC  PROBLEM 


In  general,  the  purpose  of  photogrammetry  may  be  defined  as, the  de- 
termination of  spatial  coordinates  of  objects  recorded  on  several  photo- 
graphs taken  from  different  points.  Consequently,  we  have  the  problems 
-of  the  geometric  conditions  which/  exist  for  each  individual  photogram- 
metric  camera  as  well  as • the  geometric  • stations  between  such  cameras,. 

■ The  projective  geometrical  conditions  existing  between  pairs  and' triplets 
of  photographs  taken  from  different  points  are  the  fundamental  relations 
for  the  process  of  triangulation.  The  nature  of  this  triangulation  step 
is  the  distinguishing  factor  in  photogrammetric  measuring  methods.  The 
division  is  made  according  to  the  manner  in  which  the  individual  photo- 
graphs are  being  combined  in  order  to  triangulate  the  single  points  pf 
the  model.  An  error  theory  concerning  the  final  errors  of  the  coordinates 
to  be  determined,  or,  in  other  words,  the  determination  of  model  deforma- 
tions, will  be  affected  by  the  selected  method.  The  manner  of  error  pro- 
pagation in  the  triangulation  procedure  depends  upon  whether  the  triangu- 
lation is  performed  by  intersection  or  stereo-photogrammetry,  If  is  also 
dependent  upon  whether  the  raw  material  is  obtained  by  ground  or  ji  airborne 
instrumentation,  Finally,  it  is  necessary  to  distinguish  between  a 
numerical  and  an  optical-mechanical  reduction.  In  case  of  the  latter,  the 
design  characteristics  of  the  evaluation  instrument  and  the  sequence  of 
operations  during  the  evaluation  procedure  will  influence  the  error  pro- 
pagation/ Consequently,  the  individual  triangulation  method  applied  to  the 
same  raw  material  may  lead  toil  different  results'  depending  on  the  reduction 
method  used. 

However  different  the  theories  of  the  triangulation  procedures  within 
the  different  fields  of  photogrammetry  may  appear,  it  is  evident  that  all 
methods  are  based  on  the  same  raw  material,  namely,  photographic  records 
obtained  from  photogrammetric  cameras.  It  seems  justified,  therefore,  to 
define  as  the  fundamental  procedure  for  all  photogrammetric  problems,  the 
orientation  of  photographs  in  such  a way  that  each  individual  photograph  . 
will  have  the  same  orientation  in  space  which  it  had  during  the  moment  of 
exposure.  In  this  connection,  it  is  immaterial  whether  the  orientation 
of  the  photogrammetric  camera  is  obtained  fran  relative  or  absolute  con- 
trol points  or  by  circle  readings-  or  other  instrumental  auxiliaries  e It 
is  further  ..insignificant  which  geometrical  reduction  procedure  is  being 
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applied  for  the  orientation  and  triangulation  procedures » The  actual 
photogrammetric  problem  is  basically  the  same  as  any  other  measuring 
problem:  that  is,  chiefly  a problem  of  precision  measurement  in  addition 

to  one  of  exact  mathematics.  Independent  of  the  triangulation  method  and 
evaluation  procedure,  the  fundamental  photogrammetric  problem  originates 
from  the  fact  that  during  the  reduction  procedure  the  ideal  orientation 
of  the  individual  camera  can  cply  be  apprentimated.  The 'final  coordinate 
determination  will?  therefore,  be  affected  by  the  propagation  of  the  indi- 
vidual c amera  errors  of  the  interior  and  exterior  orientation  during,  the 
process  of  triangulation,  in  addition  to  the  residual  errors  of  the  plate 
measurements  of  the  image  point  under  consideration.  ''s  " 


; In  the  following  chapters  a study  is  made  of  the  orientation  problem 
of  an  individual  photogrammetric  camera  and  the  influence  of  the  differ- 
ential changes  in  the  elements  of  orientation.  Both  the  elements  of  in- 
terior and  exterior  orientation  are  considered  and  formulas  are  derived 
which /;deal  with  the  different  possibilities  of  orienting  the  axes  of  a 
photogrammetric  camera  with  respect  to  a given  coordinate  system. 

THE  ORIENTATION  OF  A PHOTOGRAMMETRIC  CAMERA1 

(l)  The  principle  of  mathematical  perspective. 

- . ;/ 

i ’ 

i„  The  projective  relation  between  two  planes  is  expressed  by  the 

fractional  linear  equations:  £e.g.  (f>)3 


and 


where 


X 


y n 


•a^x  + b^y  + 

V + + 1 

V b2y  + c2 

V + bc7  + V 


(1) 


X,  Y denote  the  Cartesian  coordinates  in  one  plane 

x,  y the  corresponding  eoordinat/es  in  the  other  plane 

>,  ■'  ‘ ’ •'  • ' ' ■'  •* 

r In  equations  (1)  no  preference  is  expressed  for  either  of  the 
two  planes  under  Consideration.  Consequently,  there  must  exist  a second 
pair  of  formulas  which  express  the  reversed  solution  of  formulas  (1),  . 


1j  The  physical  photogrammetric  camera  is  idealised  by  assuming  bundles  of 
J.  rays  to  be  free  of  distortion  and  the  image  plane,  to  be  optically  flat 
, without  emulsion  shrinkage. 
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The  complete  symmetry  and  reversibility  of  the  coefficients  is  the  I 

direct  consequence  of  tho  reversibility  of  the  two  central  perspectives  i 

; . involved.I 2  “ j 

' " . Jj  '<;  ' ..  ' ' j 

I " • 9-  : ' T i|  J •’ 

[■.  li  ,s  ,,  f P.  Tham  has  made  some  application  of  the  reversibility  of  a central 

I •!  * l perspective  in  his  paper.  (12). 

: r : '■  •’  it  -* 


On  the  other  tend  tho  relations  between  the  coordinates  X,  Y and  x,  y, 
respectively,  can  be  stressed  by  the  foilwing  functions: 

X « f1  (x,  y,  (■*),  Kj  v,  XQ,  YQ,  II,  xp,  y , c) 

l;. 

r = f2  (x,  y/(*),  K , V,  XQ,  Y0,  II,  Xp,  yp,  c) 


or 


(W 


X - f3  (5,  y,  < , a,  05,  XQ,  YC;,  H,  ay  yp,  c) 
Y = f ^ (x,  y,  x,  a,  <o,  XQ,  IQ,  H,  xn,| y^,  c) 


P'  "P' 


and  the  reversed  relations: 


x * s1  (X,  Y,  (*),  K , v,  kQ,  Yq,  H,  xp,  e) 
Y So  (X,  Y,  (*),  K , v,  X , Y , II,  y , c) 

t ,/  OOP 


or 


/ 


a • 

!i 


(5) 


■fx  3 g^  (X,  Y,  X , a,  oj,  Xp,  Yq,  H,  xp,  c) 

y * (^»  X,  *■  > as  wj  XQ,  Yo,  H,  yp,  c) 

In  formulas  (l;)  and  (5)  we  denote  (See  Figures  1 and  2): 


X,Y  the  Cartesian  coordinates  in  the  plane  at  height  Z:  +Y  is  turned 
clockwise  for  90°  into  ,+X  as  seen  from  the  center  of/ projections 0. 

x,y  the-  corresponding  coordinates  in  the  other  plane  (Negative  or 
j)iap,ositive)  j as.  seen  from  the  center  of  projection^  in  the 
Negative- +£  is  turned  clockwise  for  90°  into'+x;  and  in  the  Die,*, 
positive  +y  is  turned  counterclockwise  for  90°  into  +x. 

‘ • ' • >j  ' * ’ 

(k)  swing  angle  of  the  arbitrarily  Oriented  x,  y » system  against 
the  direction  of  maximum  tilt  in  the  x,  y - plane 

)>  ■ ;t  •••  . 

K swing  angle  of  the  arbitrarily  oriented  X,  Y - system1'  agiinst 
the  direction  of  maximum  tilt  in  the  X,  Y - plane'1  (primary 
rotation) 

v tilt  angle  between  the  two  planes  (secondary  rotation) 

a tilt -angle  component  in  the  (X),  (Z)-plane  (primary  rotation) 

<a  tilt  angle  component  in  the  (Y),  (Z)-plane  tilted  for  a , 

(secondary  rotation) 


swing  angle  of  ths  arbitrarily  oriented  x,  y - syotem_against  the 
line  of  intersection  of  the  a,  co  - plane  with  the  x,  y-plane 


A « lina  of  intersection  of  the 
planes  IRC,  acts  as 

aXtS  Of,  tulOiiuii 

0. 0’s  centers  of  projection  in  the 

nth  mtfi  positions , rasp. 

F,F,e  centers  of/^colllneatlon  in 
.planes  I w t , resp. 

- % ii. 

1,  isocenters  tin  planes,,  i a 2, 
resp. 


FIGURE  3. 


XQ,  Y , H relate  the  center  of  projection  to  the  X,  Y-plane  at  height  Z 

x , y , c relate  the  center  of  projection  to  the  x,  y-plane 
P P 

•« 

xpand  yp  are  the  coordinates  of  the  principal  point  P;  c denotes  the 

principal  distance  - or  so  called  camera  constant;  xp,  yp  and  c are  re- 
ferred to  as  the  elements  of  the  interior  orientation.  The  other  para- 
meters denote  the  elements  of  the  exterior  orientation.  ;i 

Vi  A comparison  of  the  formulas  (l)  and  (2)  with  the  corresponding 
formulas  (4)  and  (3')  shows  that  we  have  nine  unknowns  in  the  latter  gncs, 

wnjj&y  “ "■  '5 ) •{(  . 

•’  v,  («)•,  XQ,  Yo,  H,  xp,  yp  and  c t;  ' ' 

or 

a,  ®.  * , XqJ  Yq,  H,  xp,  yp  and  c ' / 

■::'1  - ■-  ii 

. . .*  ••  jf ' . 

and  in  the  former  ones  only  eight  unknowns,  denoted  by  the  coefficients, 

V al®  a2*  bo*  \r  b2»  C1  md  c2 

t t t 1 1 1 1 1 

al’  a2J  boJ  bl’  b2*  C1  c2  y 


Consequently,  the  relations  expressed  by  formulas  (l)  and  (2)  do  not 
give  a unique  solution  but  there  are  <*>  possibilities  to  bring  the  two 
planes  in  projective  relation.  This  fact  has  its  geometrical  explanation 
in  the  set  concerning  the  rotation  qf  the  center  of  colliheatiop  [(13)  page- 
2591.  Figure  3 shows  this  principle  as  it  pertains  to  our  case.  In  crder 
to  establish  a unique  relation  between  the  two  planes,  additional  infor-. 
mation  is  needed,  v.  Gruber  in  (5)  states  that  the  unique  solution  can 
only  be  obtained  if  one  of  the  elements  of  the  interior ■ orientation  is 
given.  In  view  of  the  reversibility  of  thfe  problem  under  consideration, 
this  Statement  can  obviously  be  enlarged  ih  as  much  as  any  one  of  the 
linear  parameters  x , yj,  c,  XQ,  Y , H is  sufficient  to  make  the  solution 

unique,  in  addition  the  stilt  angle  v is  sufficient  to  fix  the  mutual 
situation  of  the  two  planes;  From  Figure  ;3  we  see  that  the  axis  of  ro- 
tation (A)  and,  consequently,  the  direction  of  maximum  tilt  in  both  planes 
is  the  same  for  all  I possible  Solutions.  Consequently,  the  angles  (/c) 
and  Kj  which  are  defined  as  the  swing  angles  of  the'  arbitrarily  oriented 
x,  y and  X,  Y - systems  against  the  direction  of  maximum . tilt  in  the 
corresponding  planes,  must  be  constant  for  all  a.  possible  solutions  and, 
therefore,  must  be  functions  of  the  coefficients  only.  v,.  Gruber  in  (5) 
gives  the,  formula  tar.  (tc)  * .a  A>  which  leads,  with  our  notations  and 

formulas  (3),  to 


cot  (■<)  = 


‘b! 


o 


a^ 


(6) 


In  photogrammetrie  problems  the  cpiera  constant,  denoted  by  c,  worild 
be  the  parameter  most  likely  to  be  given  or  at  least  obtainable  by  an  in- 
dependent camera  calibration.  ii 

■ ;7  ' 5 ■< 

lt  Halonen  has  giVen  in  (?)  two  formulas  which  express,  for  the  case 
that  c is  given,  the  H coordinate  and  the  v angle  as  functions  of  c arid 
the  coefficients.  Using  once  more  the  principle  of  reversibility  and 
formulas  (3),  we  obtain  from  Halonen' s formulas  the  modified  and  generalized 
relations*  „ 


&UU 


•2  <2 
a + b 
. coo 

,ln » ■ f ‘rT^f 
V*o  *bo 


2 V.  2 

H »o  »bo 

— « ■ I 

A r^2  TjT 
\ s w ♦ b 
V .o  o 


where  A and  A'  again  have  the  meaning  as  need  in  formulas  (3).  The  re- 
maining four  unknown  parameters  are  the  coordinates  of  the  center  of  pro- 
jection projected  into  the  X?  Y - plane  and  into  the  5,  y - plane#  denoted 
by  Xq,  Yq  and  x^s  y^,  respectively#  „ 

From  Figure  1 we  read  , 

XQ  ■ Xj  - H tan  cob  K 

: ' v (12) 

Y « X4  - H tan  * ain  K 


and  again  from  Figure  1 or  by  using  the  principle  of  reversibility 

■ if. ' _ = v . '■  !!  ■ 

i ' - xA  - c tan  ooe  (*)  'j  ;*-e. 

. : ■ a?) 

■ a ' : yD  - yi  - C tan  7 sin  (*)  ...  . r ^ 

lA  formulas  (12)  snd  (13)  X^}  Yi  and  x^,  y^  denote  the  coordinates  of  the 

isqcenters  I1  and  I in  the  X,  Y®plane  and  in  the  x,  y-plane,  respectively e 
As  ".may  be  seen  from  Figure  3 the  isocenters  remain  constant  for  all  ao  1 
possible  arrangements  and,  consequently,  the  coordinates  of  these  points 
are  functions  of  the  coefficients  only.  v.  Gruber  in  (f>)  gives  for  2^ . I . 
expressions  which  are  in  our  notation  A * 

. ao(al  * b2)  " bo(a2  * V ' „ ao(a2  “ V * bo(al  * b2) 

x. Yi r— i ‘ 

" ao  + bo  ao  + bo 


ao  +bo 


and,  consequently,  the  reversed  relations  are 


M8!  * b2-)  ° b0(*2  " bl)  - _ s0 ( a2  “ hj)  4 bp(4  * b2) 


*2  . h»S 
a0  + bo 


'2  . »2 
♦ b * 
o o 


Applying  formulas  (3)  to  (lU),  we  obtain  from  (12)  and  (13)  after  suitable 
transformations  and  using  again  the  principle  of  reversibility* 
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By  means  of  the  derived  formulas  all  elements  of  orientation  may  now  be 
determined,  namely  (at)  and  .£  are  obtained  from  formulas'  (6)  and  (7),  H 
from  formula  (9),  v from:  f ormula  (ll),  its  rotational  components  co, 

K from  formulas  (0)  and. 'the  remaining  parameters  X , Y.  x , and.y  from 
formulas  (l£)  and  (16).  ••  ° ^ P 

M ..  . ' >t 

be  spite  the  claim  of  some  authors^  that  the  prescribed  approach  give'’ 
h method  of  maximum  a c curacy  for  the  determination  of  the  elements  of  ori- 
entation, these  formulas  have  a serious  limitation.  Not  only  are  the 
formulas  correct  only  for  the  mutual  situation  . of  two  planes,  but  they  will 
fail  in  the  case  that  the  tilt  angle  v becomes  zero  or  180°.  In  practical 
photogrammetry,  however, '■ this  case  is  the "most  important  one  for  aerial 
problems.  Moreover,  it  occurs  in  special  capes  of  ground  photogrammetry 
as  well.  Even  ^though  in  such  cases  v / 0°  dr  180°,  the  tilt  angle  may 
approach  these  values,  so  that  the  numerical  computations  with  the  above 
derived  formulas  become  -unreliable.  ;j 

• . //  ti  • . . 


To  study  this  problem  further  we  will  first,'  express  the  coefficients 
as  functions  of  the  elements  of  orientation.  Omitting  the  derivation!! 
given  in  the  next  chapter  of  this  paper,  we  obtain  from  formulas  to 

(U8)  by  comparison  with  the  formulas  (l)  and  (2)  , for  the  case  of  two 
planes:  * 

? ,,  » ' II  " 

(a)  for  the  unprimed  values:  ¥ " jj  1 

_ cos (k)  tan  v tan  co  sin  K + tan  a sec  co'  cos  K 

id,  nr.ni*—.  ■ ■—  '■  ■ '■  '■  ■—«■»»  ■ ■■  1 > " ■ - ■'  ■ > — ■ ■ ■ — * 

o D D ■■ 

*’  . ■ - ;J 

, _ siA(>f)tan  v _ -tan  co  cosK+  tan  a sec  co  sin  AT 

0 TI  13  x (i?) 

a]  “ ’‘d  (-cos(tf)cos  K + sin(K)sin  X sec  v + cos  (X0 tail  v) 


IS 


n K tan  a tan  co  - cos  a:  sec  co  + (sin  Xian  01  1 cos  A tan  a sec  co)l 

J » 


e.g.  ft.  S,  Ilalonen  in  (?)  pages  lU  and  18. 


Z-ZQ  X0 

=> (-sin(#)oos  K - cos(>C)sinXaec  v + -*-7 — sin(Ar)  tan  v) 

1 D • 


Z-Z, 


-cog  >c  tan  to  tan  a - sin  /C  sec  to 

Xo  ' .1 

+ — (-cos  a:  tan  to  + sin  /c  tan  a sec  co)J 

o '« 


Z-Z  J X 

C1  “ “5 — |c(tan  v cos  K + 

L o 


X 

+ x (cos(/c)cos  K - sinCtf)  sin  K sec  v - co3(*r)tan  v) 
™ o 

+/y  (siu(<C)cos  1C  + cos(/C)sin  JC  sec  v - y-ig-  sin(#)tan  v)J 


Z-Zo  r ,>  X 

15 — jo(t«n.a-+ J5JJ-) 


* x [-sin  K tan  a tan  a + cos  fc  sec  co  - y »■ -(sin  K tan  to  + cos  K tan  a sec  oil 

•/  P U ' ’ 6-^.  J 

;‘l  ■ , ■ . ^ 

I + yp[cos  /f  tan  a tin  to  + sin  K sec  to  - ^•"^■-(-cog  K tan  to  + 3in  K tan  a sec-  «)]) 


Z-Z 


il  U Li  ^ ’A 

®2  ^ ~-g--~(-c°o (/c)sin  £ - sin(-<c)cos  H sec  v + ■£_£-  cos(A:)tan  v) 

Z-Z  - Yo  | _ COnti* 

■ -g-—  I -sin  K sec  a + yp£-  (sin  K tan  to  + cos  K tan  a sec  to)J 

Z-Z  Y 

b0  *»  — (-sin(£)sin  jf  + cos(tf)  cos  £ sec  v + sin(/c)tan  v)  a 

■ o 


Z-Z  ' Y 

" Lcos  ^ 3ec  a * z "zzr  (-C03  ^ ,t'an  w + sii1  K a sec  “M 

Z-Z  p Y 1 , 

■ "g  ■ — e(tan  v sinK  + ■ ) * 

* ■■  Y 

+ x (cop(A:)3in  K + sin(^c)cos  K sec  v cos(«)tan  v) 

P / »r  • , "o 

-■  •,  y 

+ (sin(A')sin  K - cos(/o)cos  K sec  v - g-j-  sin(/c)  tan  v)J 


, \ 
// 
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z-z 


c„  = 


2 ' IT”  1 


f 


, i 

c-(tan  co  sec  a + — ) 
L^L 


■f  £sin  K sec  a - (sin  K tan  co  + cos  K tan  ajsbb 


(1?)  cont, 
»)] 


\ - Y “ 1 

+ yp  £-cos  K sec  a - (-cos  K tan  co  * sin  /c  tan  a sec  co)]j- 


where 

D **  c - tan  v fx  cos  (<)  + y sin(&)] 


■ c - x ( tan  co  sinuK  + tan.  a sec  co  cos  K ) 

/.'  a ' P • ••  ,v» 

- / ,<r  ■ 

7 ••  % > 

+ y^  ( tan  co  cos  « - tan  a sec  co  sin  K ) 

and  (b)  for  the  primed  valxies: 

_ 1 _ cos  £ tan  v _ tan  a 


D« 


IP" 


X 


^ 1 _ sin  £ tali  v tan  co  sec  a 

O a .-Mi—  mm,m  —■—■■■■  n » ■ ■■ 

o D*7  D? 

f rt  A 

» gy  (-cos  K cos (/c)  + sin£  sin^sec  v + -2—  cos  £ tan  v) 

c *T> 

» gY  (-sec  co  cos  k + tan  a tan  'co  sin.A;:  + — tan  a) 

, x 

bx  ■ gY  (-cos  & sin(A:)sec  v -sin  £ cos  (A?)  + -E  sin  £ tan  v) 

- jj,  (-sec  a sin  k + -£  tan  co  sec  a)  ,,  .. 

r ^ , ' 

cj  ■•^7  J^,Xo(cos  K cos  (sc)  - sin  £ sin(/f)  sec  v - -E-  cos  K tan  v) 

+ Y ( sin  £ cos(tf)  + cos  K sin(/e)sec  v - -£■  sin  £ tan  v) 

u //  C 

+ ’ (Z-ZQ)  (tan  v cos  (fc)  + )j 


(18) 
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c r / xn 

c J = «•,  X (sec  co  cos  K - tan  a tan  co  sin  K - tan  a) 

1 D1  l 0 > c 

X 

+ T (sec  a sin  K - tan  co  sec  a) 

0 c 

+ (Z-&Q)(tan  a sec  co  cos  K + tan  co  sin  K + -E-Vj 

y 

a2  “ ("s*n  ^ cos(/c)sec  v - cos  K sin(*)  +-1E  CosK.  tan  v) 

c 7 

■ p,  (-tan  a tan  co  cos  k - sec  co  sin  k + -E  tan  a) 

y_ 

b2  “ C*  'cos  ^ cos(*)  sec  v “ sin  K sin(K)  + -£  sin  K tan  v) 

y 

= g-,  (sec  a cos  k + -£  tan  co  sec  a)  I 

y 

c2  3 S’1  Cx0^cos  ^ sin(ic)  + sin  K cos(/c)sec  v — 2 cos  K tan  v) 

y 

+ YQ(-cos  K cos(K)sec  v + sinK  sin(ic) E sin  K tan  v) 

y_  _ 

+ (Z-ZQ)(tan  v sin(*c)  + -E)J 

c T 7-j 

■ jj-,  XQ(tan  a tan  oo  cos  K + sec  co  sin  k - tan  a) 


+ Y (-3ec  a cos  k - -E  tan  co  sec  a) 
o'  c 


* (Z-ZQ)(-tan  co  cos  k + tan  a sec  co  sin 


>«*¥] 


where 


D*  » (Z-ZQ)  - XQ  tan  v cos  K - YQ  tan  v sin  K 


(Z-Z„)  - X„  tan  a - Y„  tan  co  sec  a 
'oo  o 


For  v approaching  zero  or  180  , or  in  other  words,  a approaching  zero 
or  180°  and  co  approaching  zero,  we  see  that  aQ,  bQ,  a£  and  approach 

zero.  Formulas  (6)  and  (7)  render  undetermined  expressions  for  (k)  and  K 

of  the  form  which  is  obviously  explained  try  the  fact  that  for  v - zero 
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i 
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or  180°  there  is  no  longer  a defined  direction  of  maximum  tilt  and,  conse- 
quently, the  swing  angles  (k)  and  K become  meaningless.  In  the  same  way 
the  formulas  (9),  (10),  (11),  (15)  and  (16)  fail  to  give  a determinate 
answer.  These  facts  lead  to  two  conclusions. 

First,  it  is  necessary  to  present  the  orientation  problem  in  such  a 
way  that  the  parameters  involved  can  be  so  arranged  that  for  the  majority 
of  practical  cases  they  maintain  their  significance  independent  of  the 
geometrical  configuration.  Such  an  approach  will  not  only  help  to  sinrolify 
the  problem  but  must  be  considered  as  necessary  for  a general  analytical 
solution.  We  satisfy  this  requirement  by  relating  the  orientation  problem 
to  such  rotations  as  expressed  by  the  a,  o,  k -system.  The  relations  to 
the  v,  K , ( k ) -system,  if  desired,  are  given  by  formulas  (8),  geometrical 
conditions  permitting.  Problems  of  aerial  photogrammetry  and  special 
problems  of  ground  photogrammetry  (for  instance,  trajectory  measurements 
where  the  optical  axes  are  within  a tilt  v from  0°  to  about  1+5°)  can 
always  be  expressed  in  the  c,  » -system.  There  are  a few  cases  of  ground 
photogrammetry  so  characterized  that  for  tilt  angles  v close  to  90°  the 
optical  axis  may  point  close  to  the  direction  of  * Y-axis  ( K~  90°  of  270°) . 
In  such  cases  a becomes  ambiguous.  This  difficulty,  however,  can  easily 
be  eliminated  by  first  rotating  the  original  XYZ-system  90°  or  270°  about 
the  Z-axis. 


Second,  we  express  the  orientation  elements  as  such  functions  of  c and 
the  coefficients  that  the  formulas  do  not  render  undetermined  values  in 
case  art,  bQ,  a£  and  become  zero.  From  formulas  (17)  and  (18)  in  con- 
nection with  formulas  (3),  the  fallowing  expressions  may  be  obtained: 


^(Vl  * boV  * + bl^p  4 ci> 

: cV(a2  + b|)  + 1 ~~  (19) 

cV(aQa2  + bQb2)  + B(a^  + b2yp  + c2) 


where 


B 


and 


1 + ao*p  + Vp 

H2B'2(a^aj_  + bU^)  + B'(ajXQ  + b£YQ  '+  c«) 

H2B'2(a'2  + b»2)  + 1 
0 o' 

H2B‘2(a£a£  + b^bp  + B'(a£XQ  + b£YQ  + c£  ) 
H2B'2(a'2+  b'2)  + 1 " 


(20) 
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where 


B*  = 


TT 


1 

a'X 
o o 


~'Tn 

o o 


Further:- 


II  * 


iA 


(21) 


and 


c ■ H 


(22) 


where 


AA'  - (a1b2  - a2Vaib2  - a£b£)  = (l+aoc.»+boc£)  =*  (1  + a^  + b£c2) 


tan  a 
tan  to 

tan  k 


a«KB' 

o 

VHBf  cos  a 
o 

Vo  “ a2 

b2  “ Vo 


K 


[ (B'H)2  + ai2] 


-1/2 


(23) 

(2U) 

(25) 


For  XQ,  Y0  and  Xp,  yp,  respectively,  formulas  (15)  and  (16)  cannot  be 
changed  simply  because  for  a » 0°  or  180°  and  a » 0°,  not  only  aQ,  bQ,  a£ 
and  b1  become  zero,  but  a^  s b2,  a|  sb|  and  b^  a a^,  bi  s a£*  Hence, 

there  are  only  four  independent  quantities  in  formulas  (1)  and  (2),  i.e. 
the  solution  is  « 3 undetermined.  Obviously,  therefore,  the  center  of 
projection  may  have  any  spatial  position.  For  c given,  H is  determined 
as  may  be  seen  from  formula  (21)  since  B and  B*  become  1.  Formulas  (19) 
and  (20)  show  that  in  such  a case  two  further  linear  parameters  must  be 
given,  e ,g . xp  and  y , in  order  to  fix  the  position  of  the  center  of  pro- 
jection in  space.  This  fact  must  be  given  serious  consideration  in  the 
computation  of  the  spatial  resection  problem  from  vertical  photographs. 

The  computed  elements  of  orientation  will  have  physical  significance  only 
if  the  elements  of  the  interior  orientation,  denoted  by  c,  and  y , are 

independently  determined  from  a camera  calibration.  Hence,  in  such  cases 
we  deal  with  six  degrees  of  freedom  only.  Consequently,  the  relations 
expressed  in  formulas  (1)  and  (2)  are  not  independent.  Two  condition 
equations  must  exist  between  the  coefficients.  The  same  problem  exists  in 
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the  case  of  ground  photogrammetrjA  where  the  position  of  the  center  of 
projection  XQ,  V 20  is  usually  independently  determined  by  geodetic1 

measurements  in  the  system  of  the  control  points*  Hence,  again  only  six 
degrees  of  freed  cm  remain.  Again  the  relations  expressed  in  formulas  (l) 
and  (2)  are  not  independent.  The  condition  equations  for  the  unprimed 
values  are: 


aobo  + + hh  ■ 0 

2 2 2 2 2 2 

ao  * 31  + ^ “ bo  ” bl  “ ^2  ° ^ 

and  for  the  primed  values: 


(26) 


4a2  + bib2  + *lc2  " 0 

a£2  + b£2  + c-J2  - ag2  - b|2  - c£2  - 0 


(27) 


t 

These  condition  equations  must  be  included  in  any  analytical  solution, 
either  unique  with  three  points  or  overdetermined  wd,th  n-points  (n  > 3). 

The  unknown  parameters  mpy  be. expressed  from  formulas  (19)  - (25) • In  the 
case  of  ground  photogrammetry  we  may  arrange  the  reference  coordinate 
system  so  that  - J - 0 aad  H - 1.  Thus  ve  obt^n  the  foll<»ing  rela- 
txpns: 


a0«i  ♦ b0bi  * °1 

a'2  + b'2  t 1 
o o 

aoa2  + b~b^  + c*> 


y * 


o 2 


+ bl2  + 1 

L + ao*p 

(al^)2  “ ^1? 


i/U 


- ejb^Kl  ♦ a^  + bjp)3 

. VU 


(a{bj  - ajbp-(l  + a0^  * 
11  * Vl  * bo°2) 


tan  "a 


tan  co  ■ b' (a'2  + l)"^4* 
o o 


U Compare  (11) . 


(28) 


(29) 


(30) 


(31) 

(32) 
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£r>' 


tan  k 


a’c' 

oc2 


aocl  " «i 


(33) 


Formulas  (28)  to  (33)  are  especially  suitable  for  the  calibration  of 
a photogr  ammo  trie  camera^  e.g.  by  means  of  star  photographs. 

(2)  The  principle  of  mathematical  projection. 

The  limitations  of  the  treatment  as  given  to  the  orientation 
problem  of  a phot ograrame trie  camera  in  the  preceding  paragraphs  are  ob- 
vious. The  assumption  that  the  problem  is  concerned  with  the  mutual 
situation  of  two  planes  is  in  general  justified  only  in  special  cases  of 
non- topographic  photogrammetry,  e.g.  an  ideal  case  is  present  when  stars 
are  taken  as  control  points.  The  solution  with  plate  constants  is  quite 
attractive  for  the  analytical  approach.  However,  such  a solution  loses 
much  of  its  usefulness  if,  as  generally  in  ground  and  aerial  photograrametric 
problems,  the  existing  condition  equations  must  be  introduced.  It  becomes 
almost  worthless  if  additional  parameters  of  the  solution  are  given,  e.g.  from 
instrumental  readings  in  the  form  of  dial  or  level  settings.  Further,  it 
will  be  shown  that  the  analytical  orientation  by  a rigorous  least  squares 
adjustment  is  by  no  means  more  complicated  if  the  elements  of  orientation 
are  U3ed  as  unknowns  instead  of  the  plate  constants.  On  the  contrary,  the 
solution  loses  its  complexity  and  offers  the  possibility  to  compute  in  the 
most  economical  way,  by  allowing  the  reduction  of  the  number  of  unknown 
parameters  by  simply  setting  the  corrections  a for  the  given  parameters 
to  zero.  The  computation  of  approximate  values  for  the  unknowns  of  the 
solution  is  simpler  if  quantities  are  used  for  the  unknowns  which  have 
easily  defined  physical  meanings  and  which  may  become  directly  available 
from  instrument  readings  or  calibration  procedures.  For  a rigorous  least 
squares  ad justment, approximate  values  become  necessary  in  both-  solutions 
(plate  constants  or  orientation  elements)  because  the  observational 
equations  are  always  non-linear . Considering  the  general  case, we  must  re- 
member that  the  determination  of  sill  nine  unknowns  cannot  be  made  if  the 
control  points  are  located  in  a plane,  regardless  of  the  number  of  such 
points  available  for  the  solution.  In  general,  we. may  say  that  the  bundle 
of  rays  which  is  being  reconstructed  by  an  analytical  or  geometrical  . 
process  will  only  be  close  to  the  original  bundle  within  the  limits  of  the 
space  defined  by  the  control  points  used  in  the  reconstruction.  In  case 
we  intend  to  compute  the  actual  physical  orientation  elements  of  the  . 
camera,  we  must  either  determine  at  least  one  of  the  nine  mentioned  para- 
meters independently  with  sufficient  accuracy  or  else  the  control  points 
must  be  sufficiently  spread  in  all  three  coordinates. 

We  introduce  now  a spatial  rectangular  Cartesian  coordinate 
system  - (X),  (l),  (Z)»in  such  a way  that  its  origin  is  the  center  of 
projection  and  its  orientation  is  consistent  with  an  arbitrarily  chosen ' 
XlZ-system.  Further,  we  consider  the  interior  and  exterior  orientation 
of  a photograph  simultaneously,  that  is,  we  consider  the  photograph  in 
its  spatial  position.  Consequently,  the  plate  coordinate  system  x,  y 
will  be  introduced  as  another  rectangular  spatial  system  with  its  origin 
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also  at  the  center  of  projection. '*  The  plate  is  now  introduced  as  a 
diapositive.  (Figure  1.) 

0 is  the  center  of  projection.  Its  coordinates  are  denoted  by  Xo>  V 

2 . In  0,  we  establish  two  rectangular  Cartesian  systems.  The  (X),  (Y), 
o 

(Z)-system,  denoted  as  the  station  system,  is  expressed  by  the  unit 
vectors,  i,  j,  k.  In  this  system,  a spatial  target  point  II  is  determined 
by  the  coordinates  (X),  (Y),  (Z),  Its  image  point  r is  denoted  by  the 
coordinates  u,  v,  w,  respectively.  The  other  Cartesian  system  expressed  by 

the  unit  vectors  x,  j,  1c  denotes  the  camera  system.  Its  orientation, 
relative  to  the  station  system,  is  expressed  by  three  elements  of  exterior 
orientation  denoted  by  two  position  angles  and  by  the  swing  angle  of  the 
fiducial  mark  system  (x,y) . The  two  position  angles  are  K and  v or  o 
and  co,  depending  on  the  arrangement  of  the  rotation  axes.  The  corresponding 
swing  angles  are  denoted  by  (k)  and  k respectively.  (Compare  Figure  1 and 
2).  The  length  of  the  plate  perpendicular  - the  principal  distance  or 

camera  constant-is  denoted  by  c.  x„  and  y are  the  coordinates  of  the 

P . P 

principal  point  P in  any  rectangular  plane  coordinate  system  established 
by  fiducial  marks  and  denoted  by  x and  y.  x and  y are  the  plane  plate 
coordinates  of  the  image  point  r in  the  oriented  xy-system  whose  origin 
is  in  the  principal  point  P. 

Further,  we  introduce  the  image  vector  r and  the  object  vector  R, 
which  are  expressed  as  follows: 

"r  ■ iu  + jv  + kw  * xx  + jy  + Icc  (34) 

and  R = i(X)  + j(Y)  + k(Z) 

formulas  (34)  express  the  entire  orientation  problem  of  an  individual 
photo gramme trie  camera.  For  the  practical  application  it  is  only 
necessary  to  derive  expressions  which  will  transform  the  triple  vectors 

i,  j,  k into  the  triple  vectors  x,  j,  1c,  and  vice  versa.  The  transforma- 
tion matrix  which  exists  between  the  two  Vector  triples  may  he  written  as 


(35) 
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5-  Compare  (8) 


>jr>* 

i 

I 


The  combined  transformation 
matrices  are:  in  the  v, 

K - system:  in  the  a,  © ~ system: 


or  in  the  K , v - system:  in  the  a,  © - system: 

i » +i  cos  K cos  v ♦ j sin  K cosv  -k  sin  v i * +i  cos  a 0 - k sin  a 

A *■.  A 

j - -i  sin  K ♦ j cos  K 0 j ■ -i  sin  a sin  © + j cos  © - k cos  a sin  © 

A „ f,  A 

k - +i  cos  K sin  v + j sin  K sinvy+k  cos  v k » +i  sin  a co3  © ♦ J sin  © +■  k cos  a cos  © 


The  reversed  relations  are;  (37) 

i » +i  cos  K cos  v - j sinK+  k cos  Ksin  v i » -4  cos  a - $ sin  a sin  © ♦ k sin  a cos  © 

j » +i  sin  K cos  v + j cosK  + k sin  Ksin  v j ■ 0 ♦ j cos  © ♦ k sin  © 

A A AAA 

k ■ -i  sin  v 0 + k cos  v k " -i  sin  a - j cos  a sin  © + k cos  a cos  © 

* 

With  formulas  (3U),  we  obtain,  in  the 
K , v - system: 

ix  « + ix  cosK  cos  v + jx  sin K cos  v - kx  sin  v 

fy  * - iy  sin  K + jy  cos  K 0 

| kc  » * ic  cos  K sin  v + jc  sinK  sin  v ♦ kc  cos  v 

i 

> 

j g,  © - system: 

| ix  * + ix  cos  a 0 - kx  sin  a 

j jy  - - iy  sin  a sin  © + jy  cos  © - ky  cos  a sin  © ;Vj 

! A , 

| kc  ■ + ic  sin  a cos  © + jc  sin  © + kc  cos  a cos  © 

i 

I 

I 
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\ 

\ 

i 

I 
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Consequently,  we  have 

r » i(x  cos K cos  \ - y sinK+  c cosK  sin  v) 

+j(x  sinKcos  v + y cosK  + c sinK  sin  v +•  k(-x  sin  v + c cos  v) 

(38) 

ocr 

r" ■*  i(x  cos  a - y sin  r sin  © + c sin  a cos  ©)  ♦ j(y  cos  © + c sin  ©) 
+k(-x  sin  a - y cos  a sin  © + c cos  a co3  ©) 

Again  from  formula  ( 3^) 
for  theK,  v - system: 

iu  ■ +iu  cosK  cos  v - ju  sinK  + ku  cosK  sin  v 
jv  ■ +iv  sinKcos  y ♦ cosK  + kv  sinK  sin  v 


kw  « -iw  sin  v 


A 

*■  kw  cos  v 


and  for  the  a,  © - system: 
iu 
jv 


iu  » +iu  cos  a - ju  sin  « sin  © ♦ ^ku  sin  a cos  © 


A 

jv  cos  © 


A 

♦ kv  sin  © 


A a a 

kw  ■ -iw  sin  a - jw  cos  a sin  © ♦ kw  cos  a cos  © 


or 


"p" a i(u  cosK  cos  V ♦ v sinKcos  v - w sin  v)  + j(-u  ainK+  v cos  K ) 
+ k(u  cosK  sin  v + v sinK  sin  v ♦ w cos  v) 


or 


(39) 

T ■ i(u  cos  a - w sin  a)  * j(-u  sin  c sin  © ♦ v cos  © - w cos  a sin  ©) 

A \ 

+ k(u  sin  a cos  w ♦ v sin  © +■  w co3  a cos  ©; 

Consequently,  again  with  formulas  (3k)  from  formulas  (38)  and  (39)  in  the 
K , v - system: 


u “ x cosK  cos  v - y sinK+  c cosK  sin  v 
v “ x sinKcos  v ♦ y cosK  + c sinK  sin  v 
w « -x  sin  v + c cos  v 


m 
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a,  co  - system: 

u = x cos  a - y sin  fl  sin  co  + c sin  a cos  co 

(liO)  cont. 

v * y cos  co  + c sin  CO 

w - -x  sin  a - y cos  a sin  co  + c cos  a cos  co 
or  correspondingly  in  the  K > v - syaoem. 

X . +u  cosK  cos  V + V sinK  cos  v - w sin  v 
y * -u  sinK  + v cos  K 

c » +u  cos  K sin  v + v sin  K sin  v + w cos  v 

* . (ia) 

and  in  the  a,  co  - system, 
x 33  cos  ex  *“  w sin  ci 

y « -u  sin  a sin  co  + v cos  co  - w cos  a sin  co 

c = +u  sin  u cos  co  + v sin  co  + w cos  a cos  co 

The  coordinates  (X)  and  (Y)  for  any  point  along  the  vector  R for  a certain 

elevation  (Z)  are: 


(X)  ■ 5 (Z) 


(Y)  - (Z) 


With  formulas  (UO)  we  obtain: 

/ \ (Z)  (x  cos  K cos  v - y sin  K + c cos  K sin  v_)_ 

(X)  = ' -x  sin  v + c cos  v 

(2)  (x  cos  a - y sin  a sin  co  ♦ c sin  a cos  co) 
— —x  sin  c '*■  y cos  fl  sin  co  + c cos  c cos  co 

( \,- % (z)  (x  sinKcos  v + v cos  K * c sin  K sin 

\*/  s " - x sin  v + c cos  v 

(Z)  (y  cos  co  + c sin  co)  _ 

“ rx  gin  a"-  y cos  fl  sin  co  + c cos  c.  cos  co 

Further,  the  object  vector  R may  be  expressed  by 


a and  consequently  (X)  * pu 

(Y)  - pv 

(Z)  • pw 


where  p is  a scale  factor 
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-»»  A — A A 

From  Figure  1 we  obtain  r • k ■ c,  and  with  R • k ■ p.  • r • k and  con* 

]J 

sequently  ji,  » — , from  formulas  (3U)  and  (37)  we  Obtain 
” c 

(X)  cosKsin  v + (7)  sin  K sin  v » (Z)  cos  v 
^ c 

(X)  sin  a cos  tt  » (T)  sin  tt  + (Z)  cos  a cos  co 


and  from  formulae  (Ul)  and  (1:3)  we  obtain  for  the  K , v - system: 

x a c(  + (X)  cos  K cos  v » (7)  sin  K cos  v - (Z)  ain  v) 
(X)  cos  K sin  v + (7)  sinK  sin  v + (Z)  cos  v 

y . c ( "(X)  sia  K + (7)  cos  K ) 

(X)  cos  K sin  v ♦ (7)  sin  K sin  v ♦ (Z)  cos  v 

and  for  the  a,  co  - system: 

m c ( +(X)  cos  a - (Z)  sin  a) 

(X)  sin  a cos  co  ♦ (7)  sin  co  ♦ (Z)  cos  a cos  <o 

a c(  ~ 00  ain  a ain  co  * (7)  cos  to  - (Z)  cos  a sin  co) 
(X)  sin  a cos  co  ♦ (7)  sin  co  + (Z)  cos  a cos  co 

Further,  from  Figure  1: 

(X)  - X - XQ 

(7)  - 7 - 7q 

(Z)  - Z - z0 

for  the  K , v - system: 

x - - (x  - xp)  cos  (x)  - (y  - yp)  sin  (*) 

y - - (x  - 5p)  sin  (*)  + (y  - yp)  cos  («) 


OM 


(W*b) 


for  the  a,  co  - system: 


(Wic) 


x - - (x  - Xp)  cos  ac  - (y  - yp)  sin  * 
y • - (x  - Xp)  sin  x +.  (y  - yp)  cos  a 
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Formulas  (U$a)  - (U6b)  are  in  agreement  with  the  corresponding  fomulas  (12a)  - (13b)  as  given  by  v,  Gruber 


31 


THE  RIGOROUS  LEAST  SQUARES  ADJUSTMENT 

The  numerical  solution  of  the  orientation  problem  and,  consequently,' 
the  least  squares  adjustment  for  an  overdetermined  solution  may  be  based 
either  on  the  formulas  (1)  or  (2),  (1*5)  or  (1*7)  and  (1*6)  or-  (1*8).  In 
order  to  derive  the  advantage  offered  by  formulas  which  are  most  useful 
with  respect  to  the  geometrical  configuration,  we  will  only  consider  the 
a,  ©,  K -system,  thus  eliminating  formulas  (1*5)  and  (1*7)  from  further 
consideration.  Furthermore,  we  must  consider  the  fact  that  the  presence 
of  more  than  one  observation  and,  therefore,  more  than  one  residual  in 
a single  observation  equation  calls  for  the  introduction  of  additional 
weighting  factors,  a fact  which  complicates  the  numerical  work  consider- 
ably. Consequently,  formulas  (1)  and  (1*6)  are  less  suitable  for  a 
rigorous  least  squares  adjustment  when  the  plate  coordinates  x and  y are 
measured.  Hence,  formulas  (2)  and  (1*8)  seem  to  be  the  most  suited  ex- 
pressions for  a least  squares  treatment  of  the  orientation  problem. 

There  is  no  special  need  for  a numerical  solution  based  cm  formulas  (2) 
which  express  the  mutual  situation  of  two  planes  only  and,  thus,  are  but 
a special  case  of  the  general  solution  which  is  expressed  by  formulas  (1*8). 
However,  for  the  sake  of  comparison  of  the  numerical  work,  both  solutions 
are  derived  in  the  following  chapters  • 


(a)  Least  Squares  adjustment  based  on  coefficients  as  given  in 
formulas  (2). 

The  observations  of  the  plate  coordinates  with  reference  to  a 
fiducial  mark-system  are  denoted  by  jf  and  corresponding  to  the  x 
and  y axes,  respectively.  The  observational  (residual)  errors  of  these 
observations  are  denoted  by  v and  v«:  hence. 


x = i + v 
y - J?'  + V 


(U9) 


The  observation  equations  are,  therefore: 

4X  * H1  *1  „ 

v =a‘X  + b'Y  + 1 ^ 

o o 

a£X  + b£I  + c£ 
v’  = a£X  + b£Y  + 1~  “ A 


(50) 


With  approximate  values*  for  the  coefficients,  we  have: 


■ ai°  + Aai 

a2  = a2°  +^a2 

a.  - a'0  + Aa£ 

b{  « b-j°  + Ab-[ 

b£  » b£°  + Ab£ 

K - +Abo 

°i  a ci°  + Aci 

c2  " c2°  + ^C2 

* Approximate  values  are 

always  denoted  by  (°) . 
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• w 


and  the  observation  equations  are: 

v(l  + a£°X  + b^Y)  * xAaj  + YAbJ  + Acj_.  - A&Aa£  - J?YA -Af 

•i 

v'(l  + a^°X  + b£°Y)  - X Aa£  + YAb£  + A c,J  - jf'XAa^  - f'YAb^  - Af ' 


(5&). 


where 

-Ajf  a a|°X  + b[°Y  + cj0,-  jf(l  + a^°X  + b^Y) 

% “ VI  . vi  •■  ..  , V'  jj  . 

,,  „ “Aj?»  - a£°X  + b£°Y  + c£°  - jf'(l  + a;°x  + b£°Y) 

\\  ...  - \ . | _ ••  ■ 

Both  observation  equations  have  the  weight  ' ij 

■ /I  ■ " '•  ..  ■■  "v-  “ 

P “ (1  + a*  °X  + b«°Y)2  ;V  * " . 

. .f'  * v ■ . V ..  - • . 

, ••  . t!%.  • i a -,••••  ••  • . 

These  observation  equations  lead  directly  to  the  corresponding  normal 
equations.  However,  to  reduce  the  numerical  computations  it  is  advisable 
to  eliminate  Ac.}  and  Ac!  by  forming  reduced  observation  equations.  * 

- x ,,  - , \ - ■ ■ ••  / ' 

We  have : " ' :r  , ..  . '/  • 


Ac|  - - Aa}  ■ i^-1  Ab-{  + "Ab^  + 

. >1 

•v;  Ac^  - - ■ Aa^  - QLI  AbJ  + Aai  + t^~-Abi  + \ 


n ^2  n 2 n 'o 
and  the  reduced  observation  equations  ares 

u . .ii 

p « AAa£  + ifAb.J "+  CAa^  + DAb^  - I 


with  the  weight  p (5k) 


p»  = AAa*  + BAb‘ ;+  C'Aa»  + D»Ab»  - L 

c.  £ 0 0 

..  '• 

. 1 .. 


I 


p ‘ (1  + ao°X  +:;bi°Yr 


33 


i JL*m 


ii  it 


where 


A * X 

n 


B 


0 - - fk  + 

? d - - jfr  * 

-L  = -a;  + -iM 


n 


Ql  a _ flX  + 

it  * n 

D»  - - f'Y  + 

-1,1  = -Ajf'  + (££.!■] 


(Mil 


p /p  ! 


v = p'/p’ 


9 


These  equations  lead  finally  *to  the  system  of  combined  normal  equations  «\ 


L J a . tf  * 

■ [pGD  + pC'D'j  -{pCL  f pC!iL5j  “ 0 


+ [pBC  *J 

* jpCC  + pG  * G *] 

[pDD  + pP'D']  -[pDL  + pD'L*]  - 0 
jpLL  pL*L*]\  - 0 

After  the  unknowns  Aa';  AbhfAah,  AUl;.  Aa‘  and'Ab*  are  determined, 
Ac£  and  A are  obtained  fron  formula!  (53),  the  residuals  p and  p*  from  « 

formulas  (5>lt)-  and  v and  v'  from  formulas  (55).  The  final  coefficients  are 
computed  frpm  formulas  (51) . During  the  computations  we  have  the  following . 
cheeks : fpAp]  + (W)  = °»  CP^P]  + [V'Dp'j  =0,  etc.  and  [vy]  + ty»V]  - 

TPPPI  + fpp*  P'l  - LpLl.oJ.  The  final  check  is  made  by  means  of  formulas  (2). 

-The  computed  x and  y values  must  be  in  complete  agreement  with  the  adjusted// 

observations  ( / + v)  and  ( £ ' + ▼ ’ ) . -If  the  A -/values  are  large,  iteration 

becomes  necessary  due  to  the  neglected  second  orddr  terms . 

The  mean  error  of  unit  weight  of  an  observed  plate  coordinate  is 


m 


Cw]  + [v'v 
2n  - TT — 


ft* 


(56) 


3k 


AaJ 

Abi 

Aa£ 

Ab'  J!  A a ' • 

d . !!  0 

Ab' 

/ L ” ;;  ' 

V-  -■  ,i"1 

\-.  ;■  , % 

[e**] 

+ jpAB^» 

0 / 

0 ! + JpAC] 

v.  + [pAp] 

• j,/- jpAlj  ,j 

^ ii  - * 

j » 0 

[ebb] 

1 

o il 

0 + [pBc| 

..  <i+H" 

/'  -H-V ; 

• iw  o 

[eM] 

+.[pABj..  +[pAGij 

+ [pad] 

/ _ 

-Hi 

r~  \V 

- 0 

I! 


The  mean  errors  of  the  coefficients  are  then  c csnpute d by  multiplying  m by 
the  square  root  of  the  corresponding  weight  coefficient  wlduh  may  be  obtained, 
during  the  reduction  of  the  normal  equation  - system,  by  applying  the  indirect ,/; 
solution  using,  e .g.  the  method  of  Chole sky,  E,  Anderson  or  H,  Wolf,  If  later 
the  elements  of  orientation"  are  computed  from  the  final  coefficients  with  ^ 
font\plas  (6)  - (11)  'and.,  (19)  “ (25 ) and  the  corresponding  mean  .errors  are  re™ 

“■  quired,  it  is  necessary  to  compete  each  of  these  as  a mean  error  of  a function 
ir  F of  the  unknowns,  pr  even  as  a n\ean  error  of  a function  of  functions  of  -the  * i 

unknowns.  The  necessary  weighting  factors  may  in  such  a case  also  be  obtained 
/ during  the  reduction  proceso  of  tile;  nomal  equations  by  carrying  additional  .. 

" dF  - columns.  If  the  Clements  of '"interior  orientation  are  given,  ce  gonorally 
in  ferial  photogrammetry,  o$  if  the'  center  ©f  projection  is  known,.'  as  is  in 
"general  the  case  in  gitouM.  photogrammetryj\  there  are  only  six  degrees  of 
freedom  present  and,  consequently,  in  addition  to  the  formulas  (50),  the 
^cisjiing  condition  equations  (27)  must  be  introduced  in  the  least  squares 
Adjustment.  By  means  of  the  approximate  values  and  the  Taylor  series!, 
neglecting  terms  of  second  and  higher, order , we  obtain  from  formulas  (27):  . 

||  a|°A aj_  + ja£°  A + b£°A  b|  + bj_°A  b£  + c£°A  c«,  + cj°A  pi  + X1  » 0 ...  f 

*■  > \ . . . '•  > ' '!  ■ ' I!  H 

" where  X..  = a'°a'°  .«•  b'°b'°  f c‘°c'°  S , : " 

" • 1 1 4V’  1 2 i 1 2 V • -•  - (57) 

and  aj°A  a£  + bj°Ab|  + c|°A  - a£°Aa£  - b£°A  b£  U^°Acj  + X?  - 0 
I 6 o o • o • o o I n 


where  X. 


■s  ;ui 


For  ease  of  computing,  the  approximate  Values  of  the  coefficients  should" 
be  chosen  in  such  a'  way  that  they  satisfy  the1  "condition  equations  (27:),  thus 
making'  both  X^  and  Xg  equal  to  zero.  This  case  only  will  be  considered. 


The  following -expressions  are  helpful  for  this  step. 


+ h * - (•  • 

D1  c2 


-bib?' 


After  suitable  transformations  we  obtain  from  formulas  (57) 

" * " - 

Aa|  - a'Ab|  + b'  Ab^+  c'  Acj  + d'  A c£  * 


A a£  ® -b » A b|  + a ' A - d ’ A c|  * c ' A c£ 


' V -<4V  + a2%20) 

a*  * — ar- s — 


ai  + 4 


, YY  14  Y 


Yfi!  I 4?4°> 

,,P  + a,°t 

4 4 i (5?) 

4' 

al  °2  " ci 


4°  f 4° 


2 * 

/'i-i  n \ , 


4°  * 4° 


5»  ■ .A 


Substituting  formulas  (58)  in  (52),  the  Ibbservaty^fi 'equations  are:  ’ a .. 

" , v(l  + a£°X  + b^°Y)  - (Y  +Xa»)Ab{  *Xb»Ab£  +'  (1  + 3fcp Ac£^‘Xd» Ac£  - jfrA^  " 

" . ..  ” > ' -jfrAb*  - > / 

1 v * (l  + &'°X  +$'0°t).  * “Xb 1 A b£  + ( Y >" Xa *.) Ab|  - Xd'AcJ  + (l  + Xc$Ac£  - jfacAa* 

% " ''-’4  . I 

: - jf'YAb*  - Ajf» 


a t['  These  observation  equations  may  be  written  as  ,,  ' T' 

| ’ ' Y ”■•  ■ ■ ".\\P  • V ...  ■■  .. 

;;  p d A A bj  + B A b£  + G A c|  + D A c£  + E A .aJ  + F A b ^ - A / 


pi  « -B  A b * > A Ah'  - B Ac'  + G Ac'  + E'  A al  + F'  A b«  - A,./' 


where 


A - 

Y + Xa« 

• 1! 

. ■ - - 

B ■ 

Xb» 

•V 

■■ 

II 

c,= 

1 + Xc» 

n 

ij 

ii 

D » 

Xd* 

<■ 

E » 

\\ 

E»  - 

X 

1 

F a 

- ft 

F'  - 

- H 

-Ajf 

= ,a|°X  + 

b.f°Y  + c«°  - 

:r 

l (1  * 

a«4  + 
0 

• K°*) 11 

-A/‘ 

. 0.. 
- <4  i 

*-b^°r  + <4° 

- f'  (l 

+ a,0X 

O 

* k**) 

•V 


i Both  observation  equations  have  the  weight 


'(1  ♦ a*°X  * b I°Y)2 


• V ;.=  ; p " . ■"  f ' ' / 1 

\ , r-  ■■■■  ' ■ " • '"^61)  " -’j  , 

▼*  “ P V P ■ ll 

The  observation  equations  (60)  lead  directly  to  the  normal  equation  system. 

Ab'  Abl  . Ac.'  4'  A cl  - Aa'  '*  jAb'  ' A/ 

1 *■  fj  1 & j o ,o  «... 

[pAA+pBB]  0 j:  [pAG+pBD]  +[pAD-pBC]  + [pA^-pBE*]  + [pAF-p3F«]  - [pA  A /-pEA.?'  ] 

, ' "m  ij  if  _ ” "II 

[pAA+pBB]’  - [pAD-pBc]  + [pAC+pBDJ  + [pBE+pAE.j]  + [pBF+pAF^]  -[pBA/+pA  A/!  ] 

./  ’ //  a 

' ..  . • . [pCC+pBD]  0 + [pCE-pDs]  11  + [pCF  -pD?.],  - [pC  A J?-pD  A /'  ] 

, / [pCC+pEo]  + [pDE-tpCE]  + [pDi’+pCF«|  .-[pDA/+pC Aj?»] 

■’  ij  ' jjpEE+pE’E*]  + JpEF+pB'F^j  -jpEAjf+pE*A  jf*j 

(62)  ..  jj  ' [pFF+pF •?]]  -[pFA jf+pF • A 


This  system  offers  an  advantage  for  numerical  computations  because  only 
twenty-two  of  the  twenty-eight  coefficients  must  be  computed.  After  the  unknowns 
Ab|,  Abg , Ac|,  Ac4|j  AaJ,  Ab^  are  determined,  Aa|  and  A are  computed  from 

/.formulas  (E>8),  and  the  final  coefficients  from  formulas  (£l) . As  a first  check  the 
computed  coefficients  are' introduced  in  the  condition  equations  (27).  In  case  the 
A -values  are  large,  iteration  becomes  necessary  before  the  condition  equations  are 
sufficiently  satisfied,  due  to  the  neglected  second  order  terms.  The  residuals 
p and  p'  are  computed  from  formulas  (60)  and  v and  v*  from  (6l).  The  adjusted 
observations  x and  y are  obtained  from  formulas  (U?)  and  tho  final  check  is  carried 
out  by  computing "a  fset  of  x and  y coordinates  with  the  final  coefficients  and 
•formulas  (2).  Thele  values  must  agree  completely  with  the  adjusted  observations 
x » (/  + v)  and  f ( /•  + y').  During  the  computing  a further  check  is  made  by 
means  of  „ 

I '•  'I 

w + [v,vii  ■ [ppp] + rpp,p ']  * [pAfAf.6] 

' and  jpApj  * jpBp’j  ’■  0,  jpDpJ  + J^pAp = 0,  etc. 


U 


The  moan  error  of  an  observed  coordinate  of  unit  weight  is 

;;  " , , . |K.  . -■  , 

, ' a , -r  \ 

...  . „ ;i  V 2n  - 6 \ " J - v - - ' 

■'  - » ■!-  * _ fl  • 

? ■ The  computation  of  the  mean  error  of  the  coefficients  must  be  carried 
out  by  "the  procedure  given  on  page  34,  ;The  elements  of  orientation  are 
computed  from  the  final  coefficients  w^ifh  formulas  (19),  (21),  (23)  - * 
(2$)  or  . (28)  - (33),  respectively . The  jmoan  errors  of  those  quantities 
must  again  be  computed  as  mean  errors  of. functions  of  the  unknowns  as  . 
explained  above'.  ' \ . 5 

. __  \ J 

(b)„  Least  Squares  adjustment  of  the  general  solution  based  on 
formulas  (U8) . „ .1  ;; 

' ' . \\  . ss  , _ h t>  ■ *; 

The  observations  of  the  plate  coordinates  with  reference  to  an 
arbitrarily  oriented  fiducial  mark  - system  are  again  denoted  by  jf  and  jf * 
corresponding  to  the  x and  y axes,  respectively.  The  observational 
s errors  of  these  observations  are  v and  v* . Hence  we  have  again  formulas' 

ii  ■ '/'  ..  " : / ' / • •"  ./  / ■■ 

ii  ;X  * / + V " ■’  fu  " ..  ' . | / 1 

f ii/  * - • ; ; ■ , I 

y - H'  + v*  .. / ■" ' 

.»■  ' •' v . • ' / . • , " , ■ /• 

From  formulas  (U8)  .we  obtain  the  observation  equations 


c [(X-X0)A  * (X-X0)B  4 (2-Z0)c]  _ 

*,+  V “ “ (X-X  )D  + (Y-I  )E  + (Z-Z  )F~~  * xr 


••  c f(X-XjA»  * (Y~YjB»  + (2-ZjC»]  _ 

+ V'  “ (X-X0)D  V (Y-Y05b  ~ (Z-Z0)F  + 7p  " G’ 


A a i'-cos  a cos  k + sin  a sin  co  sin  a:  A1  ■ -cos  a sin  k - sin  a sin  co  cos  tc 

,.jy  w // 


B »!!  4cos  co  sin  k 


B*  ■»  COO  05  COS  K 


J. 

f 

" ’!  - § 


C'  ■ sin  a cos  K + cos  a sin  co  sin  tc  C*  = sin  a sin  K - cos  a sin  o>  cos  A: 

•'.I  " ;{ 

D * sin  a cos.  to  .f 


E * sin  co 


F .«  cos  a cos  to 


l V;  if  l 


4 


\\» 


,i.  >V 


Froin  tli6  Taylor  expansion  for  the  right  hand  side  of  tho  ahovo  aquations, 
"lie  have,  ncsiootina  terras  of  second  and  hi'-hor -order:  • 

- $1  Aft  + iS’Aiio  + tcA/C+;$  Ac  * «ar  AXn  * ^ AY  + %TA^6  + 71  Ax  .-A/ 
\3a  . •<•*>  - «"*.<  -c„,  o 0 o o • lit  ox,  * ,M 

» ■ : •■'  v'  :•  ■•"■  -■  ...  .■  ■ >■  a ■ ■ :v  ,-■•  li  P 


Aa  X >S*  Af  + ^ Ac  + k0A:‘:o'  + DY0AYo  + k0A2'o  '+  v.»  AyP 

“ ti  - a •...  P 


where  - A/  « f - / -A/*  - /’°  - jf« 


' \ 


■l  « . 


and  /Ware  computed  with  formulas  (6U)  and  the  approximate  values,-  for 
■ which  ' i f . ■ v . -• " " . 


a = a0  + A a 

O A 

co  =*  co  + Aw 


yc  a >e°  + A<^ 


X = X°  + A X 
0 0 0 

Y = H°  + AY„ 

O m 0 0 


z = Z + AZn 

o o 0 


c°  + Ac 

-o 


V vAX?:  ' (66> 


•>  u*0  A “ 

y =y  + Ay 
-p  rp  Jp 


„ Making  use  of  formulas  (hU)  and  (U8),  the  observation  equations  (65)  may  be 
“ * writtlp,,  introducing  the  'approximate  values  (66)  : . ■ \i  " ' 


w 
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qQ  _ _0 

■Ou.O  / rf.o  -0\n0  _ X 1*2 


v ®»  c°B'°  - ( /' 0 - y°)£°  “ 
L ir  1 


cos  CO 


>°  x ° A a 


■ :l 


+ ( c°  sin  A?0  - 


/°-*°  ' 


£ y°)  A <K 


CO 


( /'°  - yp  a * 


-o 

- V 

* Ac 


/ 


// 

// 


c 

" /,* 


..  . / U‘ 

. H 

-c°A°  ♦ ( f .-  S?  )D°  , ( 

„ P A?X 


(67) 


0„0  / /?0  -0,  0 
-C,  B + ( l ~ X A 

_J2 AY. 


*3?  v. 

• ' 


i\  % ' 


and  V 


Cl.,0  >>  , » o -o<!„o  A 
“C  il  + (f  - XujE  - *— 


r-y! 


cos 


0 o] 
CO  X 


A a i 


/> . 

ir  vo  oos  a.  t 


jf»°  - y° 

O . * J,T>  __0\  . 


' $ 


c° 


— y jaw 


A< 


o -O  !! 


f.O  -i 

r /'  - y, 


^-Aci 


I!  . ,.C 

-c°K'p  + (j?'°  - y®)D°  • 

: *—  AI„ 

-cfeiS 

+ E AT 

„■/  \ f-  i 0 


ii  •■ 


I 


-^C.°  + (/«°-y>0  V ; 

-L- * E — azJ 

v 0 


r v 


* Ay, 

■-  A/- 


(6?)  cont. 


/'  " 


where 


and 


• \ 


X°“  -(jf°  - x°)  cos  K ° - (Hx°  - ypsin  K.  °j  - AJ?  ■*  J?°  - / 

yV-ff’  - x°)  sin  AT  ® + (K  - yp)cos  * o.  - A ■ jf* 0 - /« 

* •*  [<M>°  * <T-7>°  * (z-zo>c°] ' . -O'.  0°.  1 - 

xp  q’  *p': 


9. 


(X-X°)D°  + (Y-Y°)Bj°  + (Z-Z°)F° 


0' 


/* 


c°  [(X-Xg)A'0  ♦ (Y-Yg)B»°  + (Z.>Z°)C«0],  _Q  _ ^ + _0 
(X-X°)D°V  (Y-Y°)E°  + (Z-ZJ)F°  ^ q 7p 

O V Vv 


and' 


A ■ -cos  a cos  K + sin  a sin  co  sin  K 

■ (i  . 

A*0  * -cos  a°  sin  K 0 - sin  a0,  sin  co°  cos  K 0 


\ D°  -cos  co°  sin  K 0 


ii  ■ . >■ 


B’°  •*  •►cos  co°  cos  K° 


UO 


k 


m 


3c°  » sin  c°  cos  # 0 + cos  o°  sin  ©°  sin  X ° 

l| 

|c'°  • sin  c°  sinX':'0  -t  cos  o°  sin  to0  cos  & 0 


' Lo 


ip  * sin  o cos  & 


0 !l 


i_0  . 0 

;E  “ sin© 


The  observation  equations  may  now  be  written 


y - Ufl  + bA  © + CAJT+  d A c + e A XQ  •>  f A XQ  + g A ZQ  + AXp  “ A % 

■!'"l  " " " ■>,  ■"  ■ >1  ' J - 

v/‘  ^ a'  A # + b*  A © + c'A^  + d1  A c + e'  A XQ  + f ' A YQ  + g'  A 20'  + A jrp  - **, 


/Where  the  meaning  of  the  coefficients  a-g  'and.  a'-g1  is  given |by  the  formulas 
"f  (6f).  The . simplicity  and  symmetry  of  these  coefficients  contribute  considerably 
/ to-  the  economy  of  the  computations. 

Before  we  form  normal  equation^  we  eliminate  the  unknowns  A£p  and  A.y 


.■/ 


1 g"  _ Ml-  Cbj[  . CcJ . r,  Ldl  . Iel>  y rf]  A y Cgl . ? ..  La  ' 

axp  TAa;  tTa<0  “ irA*  - jrr Ac  - irAX0  - irAYo  ■ n AZo + 


“Aa  - ~~  A co  - Ma^-iMac  - ~^AX  - ^AY  - Waz  + 
n n n ^ n non  o.,  •<  n o n 


The  reduced  observation  equations  are  now  \ 

' v . IS  ■■  ' \ 

v » (a)  a.  a + (b)  a “ + (c)a*  + (d)  a c + (e)  A XQ  + (f ) A tQ  + (g)  A ZQ  - L 

' " ' ■'  . ■ . V ’ (70) 

v'  = (a);»A*  * (b)'Affl  + (c) 'A It.  + (d)'Ac  + (e)'AXQ  + (f)'AT0  + (g)'AZQ  - V 
where  11  ■ !!  *'  "■ 
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After  the  unknowns  A a,  A'<of  AK  , Ac,  A3L,  Ainfi;AZn  are  computed, 

we  obtain  and  Ay  ' from  formula  (69).  The  residuals  v and  V are  ,| 

computed  from  formula  a (70).  The  final  orientation  ailments  are  obtained 
from  formulas  (66)  and  the  adjusted  observations  from  formulas  (1*9)  • 

During  the  computations  we  have  the  usual  checks*  [(a)v]  [(a) 1 v*J.  * 0 , 

[jb)v]*  [(b) ’vH  s 0 > etc.,  and  (yv]  + [v ' v H - ILL.T) . The  final  check 
is  obtained  with  the  "1  Inal  orientation  elements  and,  formulas  j (6U)  which 
must  give  for  Q and  OlMhe  adjusted  observations  ♦ v)  and  (/'  ♦ v'). 

If  the  corrections  A are  large,  iteration  becomes  necessary  due  to  the 
neglected  Second  order  terms.  The  mean  error  of  an  observed  coordinate  is* 


\ 

m * + 


GnJ  ♦.  Gr'tl 

' 2n-  N “ 


(72) 


where  N is  the  number  of  unknowns  carried.  The  mean  errors  of  the  elements 
of  orientation  are  obtained  by  multiplying  m by  the  square  root  of  the 
corresponding  weighting- factors.!!  These  are  directly  obtained  from  the 
normal  equations  if  for  their  reduction  the  method  of  the  indirect  solution 
is  applied.  \ 

|i  ..  w, 

Id  general,  the  electronic  computing  devices  are  very  useful  for  any 
solution  based  on  iteration.  Consequently,  such  equipment  seems  especially 
suitable  for  the  above  derived  analytical  solution  of  the  orientation 
problem.  On  the  other  hand,  electronic  computers  are  less  Buited  to 
handle ! the  transformation  of  angles  with  their  trigonometric  functions  and 
vice  versa.  This  difficulty  can  be  solved  for  our  problem  with  the  following 
substitutions.  In  the  formulas  (67)  we  introduce  the  auxiliary  unknowns*  ! 


r - 3in  a 
s • sin  to 
t - sin  * 

Furthermore 

.0  c 

r'  ■ cos  a 

We  have  then 


and  consequently  r 
s 


r°,  + Ar 


s°  +A  s 


where 


0 . o 

r ■ sin  o 


sin  co 


(73) 


tu  +4t  ; 
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sin  K 


1 

A°  - -r 1 °t ' 

1 
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I ' 

B°  - -s'°t° 

cVU,.0 
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' 

ft 
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| ■' 
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and  the  observs 
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■ -y  i - 


J>_0.0 


,OAO 


.0  0 
J s»  * cos  CO 


1 - a0  arid  t’° 


cos. 


l-tu 


A'°  ■ -r,0t°  - r0s0t*° 
B'°  « s'°t»° 

C,0;*  r°t°  r*°s°t|0 


D°  - r°s.° 
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„o  .0  .0 

F m J>l  S' 


W I! ' 
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\ 
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where 


X0  - -if  - 5pv°  - (I?'0  - y“)t°  j 


- S)t°  * </'°  - 3?)f° 


-A / -,f>  - l 

-A?'-  f<°-  /■ 


Furthermore 


r * A a r* 


3 «*  A CO  S'  '' 


t = A !C  t< 


A r*  ’I*  - Aa  r ,j 

J’  I) 

As*  * - A co  3°  I 


A . o i 

A K t I 


From  here  the  solution  is  based  on  iterations  which  are  continued,  until  a 
pre-established  discrimination  constant  is  obtained.  Such  a value  may  be 

established  by  means  of  a lower  limit  for  the  expression  [Ta/]|  + [ [A if  j [ . 


THE  PYRAMID  METHOD  , J 

■ . ■■■  ' ■ . ..  » ■ ! ■■  2 ' . !i 

:j  The  determination  of  the  spatial  coordinates  of  the  nodal  point  and 
the  orientation  of  a single  photograph  thken  from,  an  airborne  camera  isl  . 
known  as  the  problem  of  resection  in  space . In  general,  the  elements  of 
interior  orientation  are  known.  Thus  thb  problem  is  a special  case  of  the 
general  solution  which  is  explained  in  the  preceding  paragraphs,  ..(Com- 
pare numerical  example  No.  3),  However*  several  authors  have  spent  av, 
considerable  amount  of  effort  on  theoretical  and  numerical  studies  con- 
„cemed  with  an  approach  based  jagi  determining  first  the  spatial  location 
of  the  nodal  point,  with  the  help  of  a pyramid  formed  by  the  rays  pro- 
ducing the  photograph.  .The  rotational  components  are  then : computed  as  a ^ 
s'econd  step.  The  common  opinion  resulting  from  these  studies  ,i£i  that  the 
pyramid  method  is  not  u3oful  if  there  are  observations  in  excess,  making, 
necessary  a least  Squares  treatment.  This  conclusion  is  biased  on  the 
assumption  that  ea!ch  two  of  the  given  control  points  must  he  combined  in 

an  Observation  equation"^,  thus  leading,  for  n-points,  to  ^ 1 obser- 

vation equations.  It  will  be  shown  that  this  assumption  is  incorrect 


'See,  e.g*,  (6)  page  8 or  (10)  page  66. 


and  that  only  (2n-3)  indepondont  observation  equations  exist . Hones,  a 
least  squares  solution  based  on' the  pyramid  method  may  well  be  advantageous 
in  cases  where  only  the  location  of  the  nodal  point  of  the  survey  camera 
has  to  be  determined. 


The  unique  solution  requires  computing  tho  lengths  of  the  sides  of  a 
pyramid  with"  a triang'ular  base  from  the  given  lengths  of  the  sides  of  the 
base . triang?Le  and  the  corresponding  position  angles  at  the  vertex  of  the  i 
pyramid.  The  treatment  of  this  problem  leads  to  an  equation  of  the  fourth: 
degree.  Generally,  it  is  possible  to  obtain  approximate  values  for  the  | J 
position  directly  from  the  photograph,  or,:i  in  c/se  of  strip  photographs,  jj 
by.  linear  extrapolation  of  two  preceding  nodal  point  positions . In  such  j 
a case,  the  problem  is  to  compute  corrections  AXQ,  A YQ  and  AZQ  to  the 

approximate  values  X°,  Y°  and  Z° . Hugershoff proposes  the  following  ■! 

> ..  •"  u 

method  where  these  corrections  are  found  directly  from  the  position 
angles' at  the  vertex  of  the  pyramid. 

tl  *!?  m 

■ The  camera  has  the  elements  of  interior  orientation  c,  x^  and  y^.  | 

Given  are  the  spatial  Cartesian  coordinates  X|,,  Y^,  of  the  suitably  j 

located  reference  points,  the  corresponding  plate,  measurements  J ^ and 

and  the  approximations  X°}  Y°  and  Z°  for  the  position  of  ' the  nodal  jj 

point.  ' ' ..  • ■ • „ ./  l! 


We  sets 


-XI+AX. 


■ y:  +-AY. 


+ AZ, 


Further,  "denoting  the  lengths  of  the  glides  of  the  pyramid  above  the  •! 
photograph jty  J and  the  corresponding  lengths  of  the  isides  of  the  pyramid 

above  the'  reference  ppints  by  L,  we  have;  " I!  jj 

I ..  I:  . 

L*  » (xj|  + AXq  - X i)?  + (Y°  + AY0  - 1 1)2  * (Z°  + AZQ  - Z±)2  (77) 


C.  + C . - 2L.C,  cos  ct 


r ^ 

\ 


where  a.  , is  the  spatial  angle  at  the  vertex  of  the  pyramid  and  D. . the 

_o 

slant  ground  distance  between  the  reference  ooints  i,  j . is  obtained 


(9) 


1 


1 

i 


from  the  given  coordinates  of  the  reference  points  by 

fiiJ  - (I3  - Zl)2  + <rj  - ri)2  * (Z3  - Zi)2 

Cos  is  computed  from  the  measured  plate  coordinates  and  the  elements 

of  interior  orientation.  Denoting  the  plate  distance  between  the  points 
i,  j we  have: 


cos  a 


where 


A 


- C2  * or,  - x/  * ill  - yj* 


- - <ij  -A>*  * <*j  -%>2 

Applying  the  Taylor  series  and  neglecting  second  and  higher  order  terms 
we  obtain  from  formulas  (77): 

A " Ef+  2»o  - + 2<^  - V *Yo  * 2«0  - zi>  “K 

. x°  - X.  X°  - X,  Z°  - Z.  (7S>) 

L - £?  + - AX  + -2 i A X + — - A Z 

1 i 70  O -o  O rO  O 


x2 


where  E?  « (X°  - X±)2  + (x£  - X±)2  + (Z°  - Z±) 

Inserting  formulas  (79)  into  formulas  (78)  we  obtain,  again  neglecting 
second  order  terms, 


rO  rO 

_ L.  cos  <j.  . L.  cos  a.. 

»o  - V(l  - "j  J.-  "*3>  ♦ «o  - V°  - - -jo  - - 


10 

. cos  a.  . 


E? 


'j 

cos  a 


M) 


'A  - v<«  - 3 ~ij)  ♦ - v°  - - ro 


f O *-0 

„ L.  cos  a..  L.  cos  or.. 

(Z»  - . tt)  * (2°  - Z)(l 

Li  3 

fO2  + fO2  _ g2 

+ -i g ii  _ E«  cos  a±.  = 0 


AX. 


AX. 


AZ. 


(80) 


Three  points  are  necessary  and  sufficient  to  form  three  linear  equations 

for  the  corrections  A X , AY  and  AZ  . In  case  the  A -values  are  large, 

O'  o o 

iteration  becomes  nec^sary  due  to  the  neglected  second  order -terms.  It 
is  possible  to  base  a least  squares  adjust:^-,  for  an  over-determined 
solution  on  formula  (80).  The  corresponding  observation  equations  are: 


* k - * (I°  ' x3)(Aij  - 

1—  i j _i 

[« 


+ P ■ Yi)(AiJ  -^U(JOo-  VN  - 
j»o  - ^ * K-  z3)(ii3  - 


AX„ 


AY. 


AZ. 


(81) 


+ I [^i2*  “ 5iJ*  Aij  " Bij]  + (aivi  + aivi  + ajvj  + ajvjJ 

M 

where  A.. 


» 0 


B 


2 . m 32 
ij  ‘i  “ ‘ aij. 

i - (4  - V ^ ■ % ■ V *i " • tfi  - V 

;*i 


ai ' (fi ' V 7f  ” & ■ V aj  ‘ ^ * V 

3 v- J 

and  where  £2,  3?^,  I^2  and  have  the  same  meaning  as  in  formula  (80). 


■nabiir  of  oqwtioM  r ■ 2»-3 


The  v's  and  v''s  are  the  corrections  on  the  measured  plate  coordinates 


l and  jf‘,  respectively , The  sun.  of  the  squares  of  the  v's  and  v*'s  « 

(Vvj  + [v*v*]  must  become  a minimum.  AXQ,  AYp  and  AZq  are ‘the  three 

unknowns  of  the  solution.  If  the  number  r of  the  observation  equations 
(82)  equals  u,  the  system  represents  the  unique  solution,  because  all  v 
and  v'  values  become  aero.  In  order  to  have  a least  squares  solution,  it 
is  necessary  that  r > u.  Following  Heimert's  line  of  thought  for  solving 
a problem  of  "conditioned  observations  with  unknowns",  we  may  express  the 
unknowns  as  functions  of  the  v's  and  v''s,  in  u of  the  r observation 
equations.  By  substituting  these  functions  into  tho  remaining  (r-u)  ob- 
servation equations,  the  unknowns  are  eliminated  and  we  obtain  finally  a 
system  of  (r-u)  condition  equations  between  the  v's  and  v' 's . Conse- 
quently, for  a least  squares  adjustment,  it  is  further  necessary  that 


n > (r-u) . The  mean  error  of  unit  weight  is  then  m 


Bud 


1 


r-u 


Further,  we  know  that  the  number  (r-u)  of  independent  condition  equations 
which  exists  between  the  observations  or  the  residuals  must  equal  the 
number  of  observations  in  excess.  Each  control  point  gives  ri3e  to  two 
observations  / and  and  so  we  have  for  n-points  2n  observations . Be- 
cause 3 points,  or  6 observations,  are  necessary  and  sufficient  for  the 
unique  solution,  the  number  of  observations  in  excess  equals  2n-6.  Thus 
we  have  the  relation  (r-u)  * (r-3)  * (2n-6) . Consequently,  the  pyramid 
method  for  n-points  leads  to  r - 2n-3  observation  equations . To  set  up 
these  equations,  we  may  combine  any  two  of  the  control  points  with  each 
of  the  other  control  points  according  to  formula  (82), 


Temporarily  leaving  the  outlined  approach,  we  will  first  consider  a 
direct  solution  of  the  least  squares  adjustment  based  upon  the  r obser- 
vation equations  given  with  formulas  (82).  According  to  Helmert,  such  a 
solution  may  be  obtained  by  determining  the  minimum  of  the  function 

[vv]+[v'v'3  - 2k1(a1v1+a|v|+a2v2+a^v^  . . +A^A XQ+B1 A A Z q-L;l ) 

- (83) 

- 2kr(  . . «2T2.^rnvn«-AvHArAX0.BrAI0KrAZ0-lr) 

By  setting  the  differential  quotients  of  the  function  (83)  for  v.  . ,v  , 

X n 

v^..v^,  AXQ>  AYq  and  AZ0  each  equal  to  zero,  we  obtain  (2n+3)  equations 

which,  together  with  the  r observation  equations  (formulas  (82)),  are 
sufficient  to  determine  the  three  unknowns-  the  r correlates  k,  and  the 
2n-corrections  v and  v' . 
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The  differentiation  for  v's  aid  v* ^respectively,  gives: 
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By  reducing  the  normal  equations  (98),  we  obtain  to  kr  from  the 

v-column.  After  tbe  elimination  of  these  correlates,  the  columns  headed 
by  D . T,  D aV  and  D 7 each  reduce  to  a term,  which,  when  divided  by  Tr , 

ie  equal,  to  the  w«i Siting  factor  of  the  Corresponding  unknown,  jflf  well  f 
denote  these  quotients  by  D*flQ,  TPffl  j and  [VyIj  referring  to  th<i  un-  Lib.'' 
knowns  V V and  ZqJ  respectively,  the  mean  errors  are  obtained  by  \ 


M 


Despite  the;  complexity  of  the  theory  of  this  least'  squares 
application,  the1  computation  can  be  systematically  arranged  for  hand 
and  electronic  computing  devices.  The  computational  wide  is  simplified 
by  the  fdct  that  the  auxiliary  quantities  a^..  ,,,y^  Xj,  X^  and  Xg,  as 

they  appear,, in  the  final  condition  equations  (92),  are  formed  from  the  „ 
coefficients  of  the  first  throe  observation  equations  only.  \S 

. ' f . 7 v / . : " - 

/ - j \ 

am  COMPUTATION  OF  APPROXIMATE  VALUES  OF  BIB  UNKNOWNS  \ 


The  process  of  a least  sqiares  adjustment  of  measurements  generally 
requires;  setting  up  observation  equations.  If  these  equations  are  not 
linear,  it  is  conventional  ijtb  linearize  them  by  the  application  of  the  J 
Taylor' series,  neglecting  second  and  higher  order  terms.  In  such  cases,  : 
the  neglected  term's  may  require  An  iteration  process,  depending  upon  the 
qiality  of  the  primary  approximate  values  of  the  unknowns  and  upon  the 
degree  of  convergence  of  thjs  solution;.  In  order  to  minimize  the  number  ‘ 
of  iterations,,  a set  of  approximate  yalues  is- computed  from  the  minimum  ii 
number  of jjmeasurements  neoejssayy  and  sufficient  to  solve  the  unique  case  I 
of  the  specific  solution.  jjThe  measurements  for  this  step  are  selected  i 
in  accordance  with-  optimum  !igeometrical  conditions.  - || 

: } •;!  ■„  ;;  !|\  - 

With  the  "availability  of  electronic  computing  device^  it  is  | 
possible  to  use  sin  alternates  procedure,1  namely,  to  carry  .ojjitl  instead  of 
preliminary  computation  of  ,'close  approximate  values,  a large  number  of 
iterations  in  the  final  computations!.  Such  a procedure  may  'be  economical 
in  tjismy  cases,  however,  ip1  our  case;  it  may  not;  lead  t,p  the  ^.esired  rer, 
Suit.  ^The  reasons  are  found  “in  the  jjgeometrical  characteristics  of  our’ 
prdblenj  causing  -the  solutiop  to  converge  slowly,  or  even  to-  delay  con- 
vergence by  occasional  divergence,  fit  is  evident  that  three  factors  may 
contribute  to  the  reluctance  to  converge.  The  first  of  these  is  concerned 


with  the  geometrical  coni&tiW  of  the  specific  cans  (presence  of  |ho 
’’critical  cylinder’1).  (See  below) «. ,,  The  second  factor  involves  the 
linearisation  of  the  observation  equations  which  has  the  effect  of 
forcing  the  rdtational  movements  frpm  thfe  factual  curves  into  tangential 
displacements*  The  third  factor  originates  frcrri  the  fact,  that  in  the 
solution  some  of  the  unknowns  are  mutually  dependent.  Specifically, 
spatial  translations  perpendicular  to  the  direction  of  the  optical  axis 
are  correlated  to  the  corresponding  components  of  rotation,  where  the 
degree  of  correlation  depends  on  the  opening  angle  of  the  bundle”  of 
rays  under  consideration.  These  effects  are  shown  by  the  computed  mean  » 
errors  of  the  unknowns  ifi  examples  No*  2 and  No*  3,  on  pages  72  jjand  83, 
respectively.  In  both  cases,  the  mean  errors  of  the  translations  in  the 
direction  of  the  optical  axis,  denoted  by  c and "Zj  are  small,  as  is  the 
mean,  error  of  the  swing  angle  .If  . Howevei!*,  the  linear  parameters 
perpendicular  to  the  direction  of  the  optical  axis,  denoted  by  x^, 

cj  • '/•-  > " ,j 

and  XQ,  ,Xo,  respectively,  as  well 'ag/ the  corresponding  rotational  move- 
ment so  and  ok  show  considerable  larger  mean  errors, 

f ..  '•  '■>  ■■  !l  II' 

\ ”1  ” . i i ..  ..  ■ 

From  these  considerations  we  may  conclude  that : # 

1 1)  It  1 s'  advisable  to’ introduce,  sufficiently  accurate  approximate 
values  in  the  final  solution  in  order  to  limit  the  number  of  iterations* 

2)  It  is  preferable  to  compute  these  approximate  values  by  first  | 
making  a unique  solution  for  the  throe  translations  and  then  computing  i 
the  three  rotations.  In  this  manner/  the  solution  i3  deprived  of . the  J 
compensating  effect  of  the  translations  and  rotations  and  will  in  each  ' 
of  these  steps  converge  faster,  " " " ’ 

In  aerial  problems,  the  element^  of  Interior  orientation  of  the 
camera  are  known  or  may  bo  determined  from  an  independent  cajnera  cali- 
bration. The  problem  now  is  to  compute  ’sufficiently  accurate  approxi- 
mate values  for  the  spatial  position  of  jthe  nodal  point  and  the  corres- 
ponding rotational  components.  With  thjfi  position  known,  the  rotational 
components  may  he  computed  by  the  general  formulas  (67-71)  by  setting  1 
the  A -values  for  all  translation  components  equal  to  zero.  The  compu-  ‘ 
tation  of,  the  sphtial  position  of  the  'nodal  point  is  the  well  known  prob- 
lem of  the  ui'iiquo  resection  ®n  space,  as  discussed  at  the  beginning  of 
the  preceding  chapter.  ''  . 

. ■■■  „■  . « !)  ■ „ c 

0 (I 

According  to  Unsterwaldqr  each/ set  of  three  reference  points 
determines  for  the  spatial  resection  problem  a critical  surface,  which, /is 
analytically  expressed  by\  a circular  cylinder  which  contains  the  refer- 
ence points, and  is  normalito  the  plane  pf ■ the  three  points  under  cons! der- 
ation. If  the  position  of  the  center  of  projection  is  situated  on  of  ill 
the  vicinity  of  the  "critical  cylinder"!  independent  of  the  approximate 
values  the  solution  will  lead  to  unreliable  or  even  impossible  compu- 
tational steps. 


See  (U) 


11 


i!  . .-i  if  |t  is  not  plannad’to  discuss,  within  the  scope  of  this  report, 

>!  'error;  theoretical  problems  in  connection  with  t-hoi  orientation  of  a photqr 
grametric  camera.  However,  it  3hohld  bo  mentioned  that  there  is  no 
substantial  reason  to  consider  the  existence  of  fhe  “critical  cylinder" 

«■  as  a jsSrious  limitation  in  the  practical  application  of  a 3ingle  vertical 
photograph  for  a spatial  resection.,  It  has'  becoine  an  almost  stereotyped 
phrase  in  ail  textbooks  that  the  spatial  resection  is  worthless  for 
practical  purposes,  due  tip  the  fact  that  in  vertical,  or  approximately., 
s vertical,  photographs  the;  "crytical  cylinder"  always  effect's  the  resalt, 
:;The  fact  is  that,  for  a vertical  photograph,  the  coefficient  determinant 
| of  the  observation  equations  or,  of  the  corresponding  normal  equations 
■ !!.  asympioticaliy  approaches  zero  for  a decreasing  opening  angle  a of  the 

'■  |f  surveying  lens,  '|For  an  tjniqua  solution  based  on  an  equilateral  triangle, 

this  determinant  equals,  for  a vertical  photograph,  I)  ■ 2c^  tan^  -2,  indi- 

if  ,,,!  | j f!  « - , u " *• 

l(  | eating  the  inflppnce  of  & on  the  reliability^  of  the  resection  solution, 

» l Op  the  other  hand,  for  a •»  60°  or  90°  the  determinant  V)  will  in  such 
. I cases  always  be  ^sufficiently  different  fromfzero,;;  Because. the  fore- 

mentioned  determinant  changes  sign  as  the  center  of  projection'  is  moved 
i.  from  the  in  si de  of  the  "critical  cylinder"  to  the;  outside,,  there  will  be  <• 
f a zone  of  unreliability  of  the  solution  in  the  vicinity  of  the  "critical 
ji  cylinder".  However,  it  is  readily  seen  that  for  vertical  or  approxi-  ^ 

| mately  vertical  phonographs,  taken  with  a 60°  or  90°  lens,  the  topography 
lj  of  tho  area  being  photographed  will,  seldom  cause  the  location  of  the 
ii  center  of  projection  to  be  located  on  nr  in  the  vibinity  of  the  "critical 
cylinder"  *8  • " ■■ 

..  ' ■ ••  ..  I \ ^ ■■  h.  " S 

In  problems  concerned  with  ppund -based  cameras,  the  location  of  the 
nodal  point  is  generally  known  dad’; the  three  elements  of  interior  | 
orientation  of  the  cameras  are  usually  given  with  sufficient  approxi- 
mation* Consequently,  the  problem  is  to  compute  sufficiently  accurate 
approximate  values  for  the  rotational  components.  This  may  be  done  by 
making  a unique  solution  using  the  general  formulas  (67-71 setting  the 
I A -values  for  all  translation  components  equal  to  zsto,  / In  case  the1’ 

;;  rotational  components  are  knewn  with  a sufficient  degree 11  of  approxi- 
mation from  previous  experiments  or  from  instrument,  dial  Readings,,  the 
general  solution  will  lead  directly,  by  a few  iterations,  to  the  desired 
result,  . „ .. .....  \ - J 


NUMERICAL , EXAMPLES  /? 
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Examples  1,  (Camera  calibration  from  star  iraagaa). 

, ..  ' • ' . . f . 'i  i' 

LEAST  SQUARES  ADJUSTMENT  BASED  ON  COEFFICIENTS  AS  GIVEN  IN  FORMULA  (2) 

..  (!  .•’  ()'  ..  - :l 

We  ■use  the  following  approximations  computed  from  formulas  (l8)i 

::  S " \ ” '■  A '}  ..  - ■ 
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The  approximation  values  satisfy  the  condition  equations  (27  )i 
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• ' i!'  .. 

SI 

[j 


’ • ” * a 

xli6  obSbi’VA'bion  Qcju&fcions  (60 ) scs  « 


t*V 


C&bj) 


( ACfi ) 


(Aog)  ii  (Aa0) 


-.1331705?  -.42976455  +1.14177746  + .04147660  nj.0010U343 

-.03547833  -.563510442  +1.18591035  + .05438770  +.01866734, 

#,06779447  -.16237390  +1.05373139  +.015718??  -.01039367, 

+ .07317216  -.33383422  +1.11013046  t\+ *03221836  -.00359605 

+.09683395  -.12905062  +1.0U25732U  +^01245468  -.00775451 

+.09459578  -.45982390  +1.15169J90  +.04437763  +.00877212 

ip..  'I  ..  0 • ••  • t J ■■ 

+ .15458576  - .54609038  +l»l8oi5284  +.05270321  +.02236921 

+.1913? 808  -.34136950  +1.11261631  +.03294559  -.00036325 

+.27778343  -.24642672'  +1.081^9510  +4)2378266  -.00428277 

+.34260969  -.52860226  +1.17438359  + .05101543  +.02564994 

1.42976455  -.13317059  -.04147660  +1.14177746  +.02884596 

+.56354442  -.03547833  -.05438770  +1.18591085  +.02791896 

+.16287390  ^06779447  -.01571897  +1.05373139  +.00527036 

+ .33383422  / ,07317218  ,-.03221836  +1.11013046  +.00214661 

+.12905o62  +.09683395  -.01245468  +1.04257324  -.00210930 

+ 4*5982390  + .09459578  /.0443776?,  +1 .15169390  + .00461436 

+.54609038  +.15458576  -.05270321  +1.18015284  +.00066034 

+.34136950  + .19132808  *.03294559  .+1 .11261631 1^.00807706 

+ .24642672  + .27778343  - .02378266  j +1 .08129510  -.01338374 

■•'+ .52860226 , + .34260969  - ,05101543  | +1 .17438359  i - .02412540 


(AbJ|)  -!Ajf 

(( 

+.00037595  -.00000225 
-.00186584  -.01000863 
+.00491055  -.00060210 

V\. 

-.00071037  +.00000268 
-.00587408  +.00000499 
+.00160830  +.00000034 
+ .00592575  -.00000235 


-.00020248 
-.00494660 
+ .01666818 
-.01039314 

-.00279056 
-.00249001 
+ .00042404 
-.00159802 


+ .00000763 

+ .00000886 

\v- 

+.OOOOl6j50 
+ .00002465 

+ .00002739 
-.00000209 
-.00000254 
+ .00000153 


+ .00084601  - .00000044 


+ .00017493 
-.00450210 
-.01545824 
-.01567749 


+ .00000837 

+ .00000724 
+.00000441 

+ .00000592 


-'•4-64 


//  11  ii 

■s  !,•  l/i- , vf 


The  normal  equations  (62)  are: 


+1.425234?  » 0 


ii  A I! 

+ ; 7219499  +3.2 


, Ab 

•'  - o 


+.0177304  -.0120319  +.000034509234 


+1. 4252397  -3.2484484  +.7219499:  -.0372606  -.0136234  -.000004575907 


*9.3451773 


\S  ^ 11  ‘ 


0 A if:  +.0359535  +.0138874  „+ ^000019262371, 

u \-.v 'j  ...  ||  ..  f •.  j| 

ii  +9.8451773  +.0229839  -.0438053  +.000064499348 
ff  I!  » I’ 

I +.0031415  +.0005459  +.000000983773 

n 11  ■ ■ , 

!|  I * +.0007713  -.000000281905 

/'  : , +.000000001543 


iThe  general  ablution  of  the  normal  equations  gives  the  following  weighting  cbef-  \ 
ficients*  ' ■■/!  1 » ..%  A p zi  . •.  f * \ ' .} 

' " ■ ’ l!  ji  ' . \ ‘I  :{  ■■■  ' i\ 

radl-  -5.2496817  +6.9326889  { +2.3688838  +1.0653658  "...  +75.9930602  +4.2702078 

”TT  1 ~’1|  — 1 '.j  ..  ^ ..  . t jf  , 

+6 .9326689fBBl- -i'28 . 3565122 ' j -8,837373 9 +0.3345404  -273.0771170  -21.3270208  / 

+2.3888838  -8.8373739fv4»  -2.8759711  +0.0352950  -84.7744432  / -5.0434180 

♦1.0653658  +0.3345404  +0. 0352950^  -0.5330968  +3.3129443  -10.7288684 


+75.9930602  -273.0771170  -84.7744432  +3 .3129443^3085 .6070599  +260.4505771 

" - ../..Ji  ' ! j . ■ /'  V,  II 

+4.2J02078  -21.3270208  -5.0434180  -10.7288684  +260.4505771^7]  -2309.4688222 


1) 

UTdafl  - - 2.2912184 

V ) II. 

^66]  - - 0.7301348 

fCp0]  - - 5,3250333  ). 

‘ ” -i  \\  ’ 

\[J?j  -155.5482408 

ti 

)J[TQ  * - 1.6958686' 

- -43.0569332 

The  solution  give3  the  following  unknowne:  l 

" ■■K  , _ .I  ' V;  --  ' • „ fijl 

- 0.00002460  Ac]  - - 0,00001222  A a]  - - C 

« : A b]  - - 0,00004229  A Cg  - + 0,0000078,1  f A b]  - + C 

" ..  ‘ j a , | " ; ■ 

From  formula  (58):  ;j  " 

« « 1 v 

a&1  - + 0.00004030  A a2  - 0.00002340 

. a * ' --  ■ ’ ■ I •-  ; . “ 1'- 

' - i ■ ...  ;;  v 65  ..  ...  - - L ;; 


■ V - : ■ 


J 


,Tha  final  plate  constants  (51)  ares 


.a.  - - 0.29008861 

i; 


b1*  - 0o079|6089 
*,  0.109li9li37 


/ 4 
/ . » 


0.08UU0878 


b„  - 4 OMjo^W 

* V*” 

. , ' ; d 

c2  - - OoOOOl6301 


a 0.3^113018 

® ..  Ij  n ; 

b I - * 0.09551371 

1)0 


Mean  error  of  unknowns s 

•1  A it  ■ 

4 0.00002748  V, 

4;  0 .000012 37  * m/[GaJ 

..  \V\ 

♦ 0.00000916  - 

♦ 0*00001224  10  1 

4 0.00002876  * »/,r bb]  '! 


\\  4 0.00000394  ■.  m/j  66 

V|  - ■ t / ||L  J 

4 0.00030000  “ 

■ i!  Ij  „ a.  — 

?+  o.coo2f,95i - »:fei 


where 

m ♦ 5. 


M « mean  error'  of  an  observation  of  unit  weight  after  the  adjustment  by  means 

' ~ ■(  f " 1 : V - '•  - ■ ' " 

of  formula  (63)*  !■ 

/ ‘ f 

The  final  plate  "constants  satisfy  the  condition  equations  (27  )•  ;| 

-■'ij  1/4  O.OOOOOOdh  // 


)) 

2)  - 0.00000019" 


J . \ I 

\ i 1! 


UThe  mean  errors  in  a^  and  ig  were  computed  as  mekn  errors  of  a function  of 


the  unknowns,  corresponding  to  formula  (58). 


66  i 

its 


..  "(r-'-'V  “ ^ . ;-Vi  ’it 

=-^  T,  y JJ 


/ ! J 


w , ■ 


<1 

if 

ii 

f" 

i!  % 

“ " 5 

1 . 

ii 

ii 

1 , : 

ii 

'ii 

il,  : 

» 

ii 

Si 

1! ' 

vt 

*V  ; 

r i 


-v'  ":j 

Computation  of  ihs>  residuals. by  formulas  (60~pl) 


1 

■ P 
hi crons 

" ■ P"t  V 

ii 

▼ ’ g 

. inlCroab'  S 

ii  .n 

L v 

l!  ' ' 

1 microns 

Jl 

y , 1 ’ 

+ 6.39 

+ 5.66 

/.*  d 

V 5-To 

• 2 

-10,92 , 

. +■  5.6? 

■«.  O Oo 

J *' 

If.  u.  75 

If 

3 

+ 2.30 

+ .0,83 

♦ 2 « 20 

■ il  0.79 

it 

+ 3.78 

* 6.03 

+j  3.38 

■ J 

« S.Uo 

.■  p 

: ' 5 

- 1.62 

..+  3.36  „ 

_ fl.SU  V 

A 3.20 

'i 

* $ * ' 

l.itit  | 

-8.93 

■■  . . "■  ' . 

-!1.2h' 

/;  ..  ft 

;-?•» ' 

7 

- 9.60  \ 

- 2,10 

- 8. Oil 

- 1.76 

* ■ 

« • ■ 

T • P ■ / 

J 

+ it.ie’ 

+ 3.5l  ” 

+ 3.70 

+ ,3.10 

>1 

■| 

ii 

..  , f? 

9 ' • 

+ 0.U3 

- I.it9 

t 0.39 

- 1.35// 

. j 

« 

10  ..  . 

i5-6? 

- 1.63 

H . 

+ it.  71 

. - 1.35 

pinni  Check 

s\-(jF+v)  y7{jfV) 

4 ..  H 

8th  decimi  place 


-i  - ¥ 

M V 

' , -I 


. II. 

>’  n 

!(•  U 


/ ’ " ‘ M " [vv]  j*  U08.6  i 

^ [pLLo6j-  U08.7 

Ifnal  orientation  elements  commuted  from  final  plate  constants  by;; 

” ' " !i  " ;i  ' ..  ,1  ■ I:  !: 

formulas  (28)  - (.33)  j , jj 

■ ''-U.  K - lU°  23'  5M\,  c - i 0.3007U3  H ' 

■ i;-  ti  . ■ ■ t1  . “■■■.  ; j|  ,, 

a - 19°  20*  51.8"  jL-  f .oi G00057  m 

ft  < ” ■ * - ••  i«  * i/  |i  _ <H 

a » -+5°  8»  58. it"  y - l 0,000279  m 


yp-  i 0,000279  Jn 


, Ii 
\\ 


-1 

Ii 

j +1 

!!■ 

II  0 

i -1 

0 


■ n, ' iJ  a 

. //  R I! 

. <v  • , | 

II  ' fie 


»if  . 

I. , t 


1 


Eysiiple  2.  ( Same  camera  calibration  as  Example  1) 

TOSI  mums  awnim  o«W  ohusA  sonmots  msto  w ramuu  <W) 


I 1 

.. 

j :$ 

f ;;f  jj 

! n . I 

1 .J  » 

i 3 

: ii  H ii 

1 \ -V  | 

:•  ii 

*A  1 

' s • 


I ' fs' 

1 fi 


. “j 

ii 

i 

■ t:li 


Given: 


1 ' 

- \ 


T0  ‘ 0 


L « 
ii'.  'hi 


ii„  il 
» iV  it 


"f 


Z„  - 0 

O 


and  'Iji  » 1 .. 

ii  ’ ' il 

| 3 

Approximation  values  of  the  unimjawnas 


w 

\\ 


k 0 *•  li*°  23» 


«°-  19°  23* 


V 1 

/:Ic 


O r*0 

00  » 5 


3* 


ii 

"r 


> u 

h ft 


,// 

h° 


:^°  - / o.2U7iOl8l 


•«  il 

b}< 


C°  - + O.3U2W013 


* 


I Computed  Reference 

\v.- 


' Coordinates 


1 

2 


- J 

| 

M 


X 

„ meter 
+ .1*0758996 
+.5314*6718 


« H 

I .. 
■ i\ 

I)  v 

il  I* 


meter 

\ 

K? 

• .0531*2108, 


3 I ♦ .1&U7008 
U f +.31660935 


5 

6 

7 

8 
9 

id 


-♦07298023 

+.0625U320 
+ .12239199  » | ♦♦0927250° 

+ .1*3609833  I +>*07995540 
+ ,51791371  ' t .137X9873  ) 

+.32375583  +.1801*5918 


+ ,$0132793  + * 325779U7 


is ' J 


: 11 


ii  a 


//  ■'  // 
//  Jt 

)>  if 

:4I 

' ’ ' \ ' 1 , 

••  a 

V « • 

• * "/A- 

• ■ ■) 

\ , • 

»/  • \ J 

* o? 

c *■ 

; • t.  '•/ 

0.300500  m < 1 

, , ' ' a ■ -A.  /-;J 

!l/-0\ 

/■X  » 

0,000000 

u •• 

a 

O,  O.B. 

cj-  ..  * 

' “ ’ • " ‘i  r • ■ 

0.000000 

• t ,‘A  -:v  ' 

// 

/ 

*'  ■ !! 

; ■ - 

•'  ll  ,'/  ,1  1 6 

« •? 

0*26309251* 

' D°-  + 

/ 

0.33055561  ; " 

0.961*77036 

• a ■ r;°-  + 

0,0891*7375,  . ' 

0. 00068661 

f°»  + 

0.93953576  V ; 

l 

from  (67) 

\j...  ■ 

Ml  A 

11  ••  * 

*1 " 

m 

ii  « a 

w 

" ..  a 

■ • '1 

+ .009351*25 

-.21*822757 

+1.061127  3l*  ' . 

-.1287601*6  | 

-.1911*6679 

+1.1111*2710  ...v 

- $ 

# +.2201*9851 

-.1103(821*8 

4 .981*06690  , : ij  'A 

/ - r ■■  • \ , 

ii  ■ ■ 

■'♦.01*00091*6 

-.0222^112 

+1,,oU978073 

^,20657583' 

+ .0579^*1*61 1 - 

+ V98628958 

(f<f-,0726f»33 

-.03690901! 

+1.09081*1*1*2  t ■ 

f +.00^827 

-OOSf'fs 

*.03961053 

‘if  ‘ J) 

♦1.12301073  ■ K. 

(+1.06270X1*3  „>  f 

+.06183391 

+ .1996261*1 

♦1.01*091*283 

-.1925 '4+5? 

+?.\8^9335 

If 

+1.111*1*0123  1 „ ; 

„ y . ,,  is 

a 

;i  • 
II  ■ 

, ‘ r .;|- 

o' 

\j  - 

1 i 

' » '' 

! 

!■ 

“ 

i.  ■ 

Measured  Plato 

I 

from  (67)  , ; ' j 

from  (6 7)  . 

Coordinate 

0 

~f4a 

a £j?'  **■  ii  ■ 

n 

Jj 

ft  0 ® 

jf  ■ +c 

//  ° 

J(  * -C 

; ;■  f 

V, 

t 

t- 1 

/'°-  /’  | 7 

11  , {j 

i 

+ .00264902 

-,070295m 

+ ,002560 

- .070772 

+,00008902 

+ .00047659  ii 

■\  \ 1;  , . • 

2 

..  H 

-.0348133? 

05176747 

. « i«  . . 

-.03492? 

ii  , Ij 

-.052237 

"<* 

+.00011363 

'j 

+ .00046953  1 ‘..l;  v.' 

3 

+ .0(3733262 

-.03370088 

+.067286 

-.03)4119. 

w 

-.006780 

+ , Qot>OU662 

1"'. 

+ .00041812  ■,  7 H 

: "•  | " " . ..  i \V 

ii 

:!u 

+ .011)5263 

-.00637507 

+.011358 

+.00009463 

\ » 

+,000L;0493 

; / ■ »■■■ 

5 

+.06341971 

+.01761868 

.063358 

+.017234 

+.00006171 

y 

+.00038468 

i;v  . ; If  , V* 

6 

• -.01999955  , 

./ 

= .04305923 

-.01016750 

-.020115 

1,010581 

+, 00011545* 

. + ,co°4i5M 

■i  - 1 

-.00085824 

-.01*3191 

,.fl  , ■ . ‘ 

-.001275 

+.00013177 

//  !•'  \f  

+.00041676  /„ 

/ N rj  ■ ' ! 

1 , ■■  '• 

8 

/ 

♦ ,0012  32-39  r 

+.02533916 

+.001122 

*.021*948 

+.00011039 

+ ,00039116  ;7"  ,7  ^ “ 

...  \ /'H-  . H ; ■ 

'■!■  9 

+.01842761 

+.05762827 

+.018325 

+.057266 

+ ,00010261/' 

..  . / 

+ .|pQ36??7/|/  :i 

....  . ‘V  jj  i- 

10 

-.05100U0U 

+ .0485.0096 

-.05116U 

+.048123 

♦.OOOI5996 

. i - 7 * , 

+.00037796 

. )/  : . » ... 

" - V 

y° 

„ a 

. "■  ' •( 

Vi  •'  • 

..  t t . . 

.7  b "/ 

>.v.  .v« 

n 

//  , d/ 

1 

from  (67) 

+ .011*89596  -.06875007 

*> 

+.29607230 

froril  (67)  and  (68) 

+.07525270  +.00881537 

■1  "XI 

II  . ’ ' ‘ 

i' 

i 

2 

+.04658162 

-.04149692 

+.29992023 

+.06983916 

11585148 

r ■ 

it 

3 

-.05685054 

-.04937051 

+ .30561617 1 

♦.0&570902 

• ' ii 

+.22406862 

4 

-.00951003 1 

-.00902617 

+ .29084489 

+.07499042 ' 

+.03811191 

..."  I „ 

5 

-.06580847 

+.00131246 

l! 

+.30217008 

\ ♦, 07436965 

+.21104719 

"! 

» 6 " 

♦,02189836 

-.00488075 

+.2922J480 

+.07432181 

-.06655424 
-.14329195  ;• 

7 

+ .04192?, 75 

+.00986492 

,♦.29597330 

" (/" 

, +.07606020 

■y 

..  \ 

8 

- .00748821 

+ .02423878 

f + .28767689 

> . .jj! 

+.07454723 

# +.00410113! 

. II  h » 11 

l ) } 

9 

-.03216534  / 

+1.05124437 

' + .28672184  ir 

”+.07150417 

+^06132 M 

if  ti  ‘j // 

if 

10 

+ .03735731 

+ .05965Q54 

+.29188918 

| +.08477116. 

r r*i697#8  j 

' ,1  !iS 

. I.  , 

- i‘:  11  - , 
!■’  ’!  - 

.■69*  V1 

i " ■, 

/•  « v ■ 

n ii 

- "•  ' ‘ - 1 

-i 

(Jt- 

l».v 

i:  ' , ; 

i . . / 

' ! -I , l '■ 

II  j 

11  :i 

II 

0 \ /.  -■ 

• ■ Ii 

-l 

\i  ii 

!l  « 

' "/ ' ' ii  i 

it- 

\ ^ P ''  "■■' 

ii  .. 

a jj 

5 /(  " 

. ' ' ii'  - , _ 

- iiii,,; 

„.  ; til  v ,,  ■ 

»-ti'  " 

- i* 

1 " r ■J7, 

;1  I - ; ./ 

'■  y ‘ ; , 

/ ' 7‘.  . 

• jj .. 

' ‘ . ir  i! 

jj  ' - 

,!  - 1 ■»  : 
' - Ii-  - ■ v' 

. 

» 

...-■..■ii 

f ; . '•* 

.sf55* , 

T'".:  nO"  <«u 

tf  f ..  \ 

,je. 

...  ...  Z 

|i  ' ' - I'd:-..  . V 

1 .gjggjji, 

• \ 

Li]  m:  + 0.0021*077$ 

n ,i 

[s3m  + 6 ,001$ 6378 
n 

LM  - + 0.00$20392 
n 

- 0*0080121*8 

''  n: 

CMS-  - 0 loOQ102$8 

n 

i [MV-  0.0001*11$$ 

n 

- • ;'i  f;  11 

The  reduced  observatiqa  equations  from  (70)  are: 

. ' !i  1 ' 

..  ffA  il  -+ 

- T • :■  n r’7 

/ \ II  /m  \ 11  / \ 


; . • ■ ; - 

(a) 

9 (b) 

’(0) 

i/  • a 

(d)  1 

-L 

■ „ 

1 

A a 

!»  A a* 

A 1 !! 

A ic 

Ae 

i ; , ..  11  “ 

Aa  “T3C 

Aw  r isb 

// 

1 ; , ■ , l 

+.115632 

-.088^85  / 

+.6788766 

II 

+ . 0036111+5 : 

ir, 00001356 

:f  ' i : ' 2 

+.5001+25 

\\  / ti 

-.629739  ' \ 

+ .1+935972 

-.1210551+0 

+.00001105 

i il..,  „ 

If  » 3 

; (j.. 

+l.ij)70019  > 

+.957217 

,+'.3129313 

+ .218861+70,1 

■ ' 11 

1 00005596 

i ■ '•*  ■ ' 1+ 

; h f: 

-.107109 

-.111+613 

+.0396732 

+.01290799,, 

Jl 

-100000795 

' '5 

■ il  ;/ 

! +.725110 

-.176690 

w 

-.200261+3 

+ .20581+337 

-.00001+087 

| ’ ..  6 

i -.26I1118 

-.1811+71+ 

+.0775575  . 

-.07175816 

♦, 00001287 

4-  7 

1 +.1057U0 

it 

-.007635 

, -.0151+951 

-.11+81+9587 

+.00002919 

•s’.  8 

■■  n w 

n 

3 -.723909 

-.158932  * 

-.2771+691 

/ -.00110279 

+.00000781 

’■  -1  9 

? -.8i9iaii  ' 

-.1+63238 

-.6003602 

H 

+.05611921+ 

♦.00000003 

o 

iH 

'-,302680  ■' 

+.8631+61 

-.5090871  | 

-.171+931+50 

/ 

+.00005738 

, ' !■■  ■'  ' i 

♦,252976 

.il 

-I.IOI53U 

+.0108521+  1 

-.22591567 

+.0000650I+ 

”1”  2 

+.36997U 

-.20811+9 

rl  1} 

-.3637715  ... 

? -.161+25863  ; 

+.00005798 

"j  ; t;-"  >'.  3 

-.1503511 

-.01+6961+ 

+.6576881+  } v 

.// 

-.101+13687 

+.00000657 

« 

-.0351*09 

+ .1+87587 

+.0988885' 

11 

-.0132021+0 

-.00000662 

m 

• if  ’ $ 

: . . ■ 5 

" ! . i;  vf;  6. 

” i . //  a 

1 "+,8|39l+9  | 

+.1+99028 

+.6185593 

+.06661+369 

-.00002687 

-.222933 

!!!:■"  " 

+.1+90209  •!! 

, • »1 

-.2156333 

-,02582279 

, , ,;+. 00000195 

7 -.1+91X28  +.50951+0  I -.1+1+62301  + .005l56ijl+  +.00000521  ! 

' | il 

8 -,bU3^0  r +, 30233V  | -.0033139  +*09233581  -.0000203 9 1 

■'  \ S _ ..  ; 

9 +.66ll+t>9  -.1+76011+  1 , +.1686383  +.19978709  ; -.00001+928 

“ ; . ' * " ■"  • i " 

10  -1.17U670  -.1+5601:2  ,,  -.5256782  . +.1691+1335  -.00003359 r 


“I1!  - f v.  J!  . • . . <’ 

fi  r~The  decimal  point  in  the  coefficients  (&)•♦* (c)jwas  moved  fcr  the 


Aa 

+6.1*991619 


+.7961320 

,+5ii55ot*k 


■•j’i  - 


means  of  (71) 

are: 

u ' . n 1 " j 

• 

Ak 

;i  Ac.. 

'•i  ■ Ij  '«  j 

“1  : 3 I 

+2.6017281 

■ ,,  1!  5 

+ .2070662 

\\  (j 

-.000086182255 

+.01*33678 

it 

^ +.21*281*37  8 ' 

(l 

-.0000698981*31* 

+3.051*7725 

it 

■ j 

-.000070911*325 

'»j 

+.338290(1 

-.00008165861? 

1 

ti 

■ l!  ■'  : ' 

11 

'■  • 1 

+ .0000000219U7 
» ““ “ ■* 

, The  veighting  co 6thcl<Ji  &°»  the  g^ril  solution  oC  the  ncml  actions 


faVV  -^9965  . +.0300362  . ' +.2582353  +.1283995 

+ .0300361  [b  ' -.201*1*813  -,0226786  +.12^0?9 

> / p , fj  if  1 ft 

+.2082353  -.0226786 |y  Yj.  -.5(^3872  -.1111798 

+.1283995  +.1281*029  -.1111798  [ee] - -3 .1266057 


. Il  *,  V,  ■ n 

The  solution  gives  the  following  unknowns: 
Aa  - - 0.00062369  - - 2*  08.65" 

\v\.  ••  

Aw  » + 0.0002827^  - +58.32" 

A k - + 0.000281*86  « +58.76"  „ 

Ac  - + 0.00021*317  » 

Fran*  formula  (69):  a 

l , '•  .. 

Ax  - + 0.00005788  Jj’ 

> { p ■ ■ * *'  I a 

Av  - - 0.00027921  I .1 


L . J 


,/[aaj®  + 49.5 
/[ftp]-  + U5.2| 

7,1 

/5i]-  + 1.8 


l!  '-jj' 


! 


The  final  orientation  elements  from 


Mean  errors  of  unknowns 


k • l!;°  23'  $8.8'* 
' a - 19°  ,20 ' flJx' 
©«  5°  08*  §8 *3" 

5 

f c - 0,300743 
| Xp»  +0  looo58 
•5.1y  * -0 .000279 


7.9  - pftvyrrYl 


t55 .1  pm  V [aal 


50.4  - WQ#J 

+ ,000009?  =>  m /M 
i ,OOOG78."1?  , 
t .000078  12 


where 

m » t £>,4n 


p ® 206265 


n|ean  errbr  of  an  observation  of  unit  weight  after  the  adjustment  by 

. L» \ . ! ■ " - ..  \ 


formula  ( 72 ) . v, 

I ■ ' ii 


Computations  oj[  the  ■ residuals  by  formula  (68)  f 


microns 


..11  + 5.69  J 

■ :v.  it  ' 

■/I  - 9.23  / 


| .+  2.21 
j + 3.40 

I r. 1J5h 

*!”- 1.24 
;f  - 8,04 

/ + 3.70 
+ 0.38 
+ 4.6? 


microns 


+ 5.72 
+ 4.78 
+ 0.7 9 

- 5.41  y 
+ .3.17 

- i.74 

+ 3.10 
-1.36 
- 1.33 


Final  Check  by  (64) 

g~(jt  + v)  o'-a’+V) 

8th  decimal  place 


’’  j1  V - 
ill " \ 


■0  8 

1 ' 11  i 

*'i  tr  !! 

' -10 


[w]  e 409.2  [pI.L.6]  - 409.. 6 


The  mean  errors 


:0r3  in  x and  y wore  computed  as  mean  errors  of  h function 

^ p . i it  ■ '■  ii  < 


of  the  unknowns  corresponding  to  formulas  (69). 


II  ,11 


■ II 


A comparison  of  the  results:! obtained  from  examples  No.  1 And  No*  2 
...  '■  '|V  ' "■  « „ ' ..  * ■ 

shows  a sufficient  agreement  between  the  two  methods . ; ■ • 


FINAL  RESULT 


Method  !,A" 


Method  "B" 


" 

based  on  formulas  (2) 

1 based  on 

formulas';  (1*8) 

u 
. 4 

a 

+19°  20  ’ 01  .Cj- 

" +19° 

| . 'ft 

20  01  Jk 

n 

CO 

+ 5°  -ba*  08Ui1 

• ■■  \ • 

V + 0° 

08  08.3' 

23*  08 .8™ ■ ■ 1 

//■  ■ »• 

K 

+ll*°  23*  08.8", 

■ 

* +H*° 

c 

.+  0.30071*3  ; ‘-’i  . 

a 

1 + 0.30071*3 

h 

+ 0 .000007 

V+  0 .000008  . ..  ! ' 

«i  - 0.000279 

ii 

ii 

■ i U 0.000279  ■. 

'.  1 M • 

..  " ' if  i|  » 

j|  V. 

SpIDUAIS 

Method  "A”’ 

II  ' 
...  ■.  t/ 

' ‘-L. 

Method  *'B«  , 

I 

\S  • . 

V V*  „ f 

* 

v 

f.'.-.r  v1 

n ; ; x\ 

microns  microns 

II 

microns 

i microns 

i 

+ 0*7  + 0.7 

ii 

+ 0.7 

j + 0.7 

l ij  ..  ■ ■- 

- 9.2  + U.8 

l! 

\ jj  n 

-9.2 

+ 1*.8 

+ 2.2  f 0.8  ., 

" • ii .. 

••  -L  if  - 

+ 2.2 

1 ’ + 0.8 

(i 

+ 3*  -"0.U 

' //  y ■ " ' 

ii  i ■■ 

+ 3.U 

! -0J* 

' - 1.0  ♦ 3.2 

• * - ii  ' 11 

il  . 

* ^ 
- 1.0 

..  ;;  •+  3.2 

ir 

- 1.2  -7.7 

- 1.2 

' ; 1 •-  7.7 

It 

a 

- 8 .0  “1.8 

— 8.0 

- 1.7 

'!  + 3.7  +3.1 

: + 3.7 

+ 3.1 

fj 

+ Oj*  „ - 1 fb 

»! 

+■■  0 .1* 

!,  - 1J* 

+ U.7 

+ u .? 

J - 1.3 

[yvj  - 1*08,6 

[wj  -1*09.2 

il 


!’  p , 

i»\V  ft’ 


>l  A 

• mV- 


n... 


rjop^aplo  3.  (A  spatial  resection  flpon  an  approximately;  vertical  photograph.) 

LEAST  SCARES- ADJUSTmT  OF  THE  OEM^aL  SOLUTION  BASED  ON,  FORMULA  (U8) 


i .. 

« 


■\\ 


, The  computation  follows  purposely  such  numerical  :jsteps  as  are  suitable 
for  electronic  computing  devices,  given  in  fprmulas  (73)  to  .(75)# 

Elements  of  interior  orientation  as  obtained  from  a,  camera  calibra- 
tion, and  assumed  as. free  of  errors:  • .....j 


q - 0.16999 


X = O.OOOOO  i „ 
P 

V « C,  0(3000  « 

p n ; 


* « 


* :l 


■"  • ,/■ 

First  approximations  of  the  unknowns  taken  directly  from  the.  photograph : 

X?  - S5tta 


1°  ■ 6tan 


Zo  ' 


k 0 ■ 66° 
a°“  180° 


0° 


'I  •••. 
.1  .'I 


v\ 


Auxiliary  unknowns  from  (73)? 

il  ••  ■...  - | 

1 r°  - sin  a0  w,0„ 


co 


'o  o - ; 

r ■ cos  a “ - I,' 


: \ 


m 


ffi  * sin  a>°  » 0 


s*  COS  C0°  I1! 


t®  ■ sin  k°  ’• 


Auxiliaries' from  formulas  (73): 
" - | 

A°  - +.U067366U  jjV  5 


.1  ll 

tV5"  cos  k°  » + » 1+0673661* 


B°  - -.91351*51*6 


■ S 


; jj  » 

::  a » 

■ II 


' ji 

Ifl  J •>  -ts  • 

1 1 if 

II 


// 

ti 

ti . 


« 


1 ‘iii 


75 


.Ji  '«• 


?0“ 

//  • • 

+ .9135H5U6. 

D°  » 0 

I ' 

•o- 

+ .1+067366 U il 

E0  0 

,0» 

0 / 

7 ..  ■ 

F°  * » 1 

W 

fl 


■ ,1-  '.J 


Coordinates  of  ground 

,f  ji  , 

1 control  points 


•from  (<$,?) 
■ n 


„ '■»  ■ Ml  . . • Jr 

1218.4$  290.8  I4JU6.I  +590.893162  *468.627028  + 2153.90  +.04663445  -,03698496 

>■  I ■ ■ ....  I " s " '*  1;  . ..  V 

750.95  125.8  435,9  +551.478803  -25.567020  +2164.10  +.0433186$  +.00200829 

T'  ' 1:  ' 3 ■ £ \ " 11  ■ V.  •*  - " * ' ,d  s? 

273.20  ! 0 UliO.l  +47^.084392  -513.180533  +2159.90  +.03715432'  +.0403887f 

■ ’ r »;  i ' 4 ;!  ,i:  •'  ||. 

1094.2$  738.4  448.4  +131.473543  . +537.220002  +2151.6C  +.01038724  -.04244378 

“ * I (j 

621.15  622.6  437.6;/  *4.i.83$003  +57.921542  +2162.40  +.00352456  -.004553!$. 

3.20.55  392.1  442.4  +511.794869  -493.152111  +2157.60  +.00408074  +.03885379 

II  fj  •!  . i«  . II 

929.95  1177.3  470.0  -33^.308389  *565.641194  +2130.00  -.02683994  -.0451424? 

446.05  1060.4  449*7  -426.334785  +76.029033  +2350.30  -.03370351  -.00601041 

0 934.8  439.0  -493.018354  -382.544042  +2161,00 ; -.03878213::  +.03009193 


Meador ed  plate 
coordinates 


t 

+ .45978 


-Ajr- 

i fc . ( 


+.036240  +.ooo65645  +.00074496  -*05275538  -.027 55955 


+.043105 

-.002307 

+.00021365 

+.00029871 

-.01578462.  ,, 

-.04039040 

+.037465 

-.040100 

f* 00031068 

-.00028873 

+.02178492 

-.05036974 

+.010086  / 

+,041152? 

* +.00030124 

+.00129178 

-.04299919 

+.00777425 

. :w,  ' i( " 

+.003815  | 

+.003612 

-.00029044 

• -rt\  • 

*,00694131 

-.00559322 

-,dbl3678$ 

■ ( ;i 

+.004931  i| 

-#639019 

-.00085026 

+,00016521 

+.03383492’ 

-.01953126 

-.02670'?  ' 

+*043294 

-.00013094 

+.00184842 

-.03032287 

It  - 

. +.04288658 

-.032878 

+.004552 

1) 

r. 00082551 

4 ' ,Y 

+.ooi4584l 

J +.00821767 

!i  *! 

+.03323434 

i-. 037327  | 

-,030990 

-.00145513 

+.00089807 

+.04326446 

+.02318975 

[All]  + .00269162.  [A  1 ’j  007  35814 

’!  ' , JiAijl  >'iai':i.  . 

discrimination  factor  jfzn 


.0QO55632 


■V  1 
I 76  «i 

r;  - >1  ■ 

i!  ;?  .. 

• u * : [! 

4 ? 


■ 11.  ;; 

' ll  f. 

'•!  ;? 

-Wl  ' I’n' . ‘ 


g 


;1  •• 


T*  " V 

J \\ 


©**f; 


5 

ll 

\\ 

i!  - h 

+*08.361388 

+ 

* 

.162851418 

-.00003210  ; 

+.00007210 

■ - l' 

-,00002165 

; ' 

\\ 

+.07316356  A 

•» 

*i 

>16558630 

-.00003195  .1 

+.00007176 

. :i  !! 

-.00002002 

+.06437968 

+.1663021)9 

' 00003201 

+.000O7190 

..  ,j.  • ii 

-.00001720 

■ .11  ‘ •*"  • ' 

■'ll 

I +.07176863 

+.15431855 

-.00003213  * 

+.00007218 

•'i 

-.00000483 

| -■  fi  ■.  ;* 

\ • 

+.l>6925713: 

+.15532195 

-.00003197  ;i 

+.00007182 

-.00000163 

a / 

+.06832893 

1: 

+I155W5 

-.00003205  “/ 

:+ .00007198 

OOOOOlEp 

// 

V//'  ‘ • •• 

+.06435344 

i 

+ .l!62064o6 

-.00003214:6 

+,00007291 

*.00001260 

.»• '1  ■■■ 

\t  / 

y; 

/ 

1:  ^ " ' ■ 

+.07077046 

+.16188288 

-.00003215 

Jj.  0pij)0?222 

+.00001567 

/ 

S 

i 

is 

" +.07901167 

! 

. . 

+.16058U18 

-.00003199 

;+.  00b07l86 

ti 

00003795 

If  i 

1 ' 

;!  •• 

/ •• 

i!  • w ■ 

it  * 

■ i 1. 

■ ift 

T1 

1 /' 

,/  Q>  if  " iu 

! ........ 

b 1 

" f ..  ' | 

«.*  Ill 

« M 

■J, 
••  % 1 

+-  ■ ■/.  + 

. ■ s +.1667716#  „ 

-.0631UU99 

-.00067210 

-.00001717 

V 

1; 

■ ' ;;i  " x. 
|;  ,<■ 

+.15510711 

069618, 3U 

-.00007176. 

+ .00000093 

si 

V 

“ ■ ip-  . 

Hi  1 

+.160U695? 

-.08110875 

400007190 

({ 

+ .0001870, 

# i 

\\ 

\ ■ fi  5 

11.  +.16602980 

07108226 

-.00(007218 

-.00501973 

! i 

■1  f\ 

) y.  i; 

]/  11 

Vj  +.15544341 

-.06910452 

-.00007182 

-.00000211 

»// ' j 

„ 1 

!i  . : +.16302707 

-.07360531 

-.00007198 

/+.  000  01801 

£ 

i 

W' 

(r 

« 

+.16 334611 

-.06052850 

-.00007291 

' // 

-.000*02119 

*■ 

[•■-»,  , ■ 7 

. +.15500303!  .. 

-.07031624 

".-.00007222 

• V " : 

-.OOOWO280 

, , * 

him  •» 

ii  . • s.- 

+,16295234 

-l\  *• 

-.06503607 

’i+t|0007l86 . 

'Hi 

+ . 0000139 1 

it  7 

w • (( 

* ws  .•? 


V “ 


• >! 


■ ' p 
tr 


r".i 

’ rf 


')  ■ 

» : a 


■Rasaiii 


r ■ ) 

ft 


v.  :j 

.//  ‘j]  ('  ■'  • - *'  !|  \ -j  ^ 

/the  observation  equatidns •( 74 )>are : 

> 

- ii  . 

n >1 

- • ‘ ' ; ' " . 

!! " ■ ‘1  v 

< 

1 

A2e 

|-  4. 

. I Ag  'v 

Ac 0 Ak  V » ^ 

w 

' -■  U , VA..;  ; - 

i 

i 4- 

• „ • :» 

• // 

AXn 
a 0 

AY 

0 

i ' ' ‘ 

ii 

. ’ ..  . . 4 A\  ‘ 

■jj  •,  :>.* 

T6m  ~ 

TSM5 

; “ lO60O 

j\  ■ ■ V,  ’ ■'// 

i 

+.08361388 

+>16285418  -.3698496 

-.3210,, 

+ .7,210 

-.2165 

+.00065645  ’j 

1 • y 

!/ 
| / 

+.07316356 

+.16558630  +.0200829 

(1 

-.3195/ 

+ .7176 

: -.2002 

'.1 

+ .00021365'  / 

+.06437968 

+.16630279  +.4038873 

-.3201 

+.7190 

-.1720 

‘ -.00031068  '{  , 

[ *1.  • 
1‘ 

+.07176863 

+.15481855  -.4244378 

-13213 

+.7218) 

-.0483 

+.00030124 

• <.« 
I-T- 
! *■'. 

! .fljy 

+;o69257i| 

+.15532195  -.0455331 

3197  ij 

+.7182 

-,0163 

-.00029044 

krv 

j ir 

+.06832893 

+.15576245!  3885379 

-.3205 

+ i'?198 

-toi89, 
f +.1260? 

, ii 

-.00085026  \ „ v 

,,+.06435.344 

" "/  v •’ 

+.,16206406  " -.4514242 

-.3246 

+ .7291 

■V  ■ 1 \ !»/ 

-.00013094 

I i 

II  0: 

+.O7Q77046 

;i + .16188288  -.0601041  ' 

-.3215 

+.722& 

+.1567'! 

’;|  . ('.fj 

-.00082551  \ 

J +.07901167 

1 +.16058418  +.3009193 

-.3199, 

+.7186 

"+.1795 

-.00145513  >Vvv  A 

• 

' ( t +.166773.65 

■'-.063144 99 1 +.4663445  > 

,-.7216 

-.3210 

-.1717 

+.00074496' 

' 

+.15510711 

-.06961834  +.4331865 

-.7176 

-.3195 

+.0093;' 

+.00029871 

+«16o46957 

, -.08110875  +.3715452 

-.7190 

-.3201 

+.1870 

-.00028873  J ■ 

*3-: 

+. 1660298b 

-.07108226  +.1038724 

-.7218 

43213 

,<1 

-.1973 

+ .00125178  . |'” 

//  & 

+.15544341 

-.06910452  +.0352456 

-.7152' 

-.3197 

-.0211 

+.00094131  " 4 

V.  " ’ - i • l'  jj  ij  ’.v  *'  H N 

+.16302707  -.07360551  +»0408,0?U  -.3205  +.1801  +.00016521 

" ••  ■ II  l "‘  \ S'  , " ; 

-x  +.16334611  -.08052850  -.2683994  -.7:291  -.3246  -.2 119  +.00184842 

■ ■ f >? ■ ••  - 4 - *■ ' * ■ u ■ 

r +.15500305  -.07031624  -.3370351  -.7222 ' -.3215  -.0290  +.00145841 

..  « ».  " j 8 j " »■  i \ 

+,16295234  -.06503607  -.3878213  -.7186  -.31 99  +.1393  +.00089807, 

..  ' I!  „ ' ■■■  ' »■  : v h ■ ..  7 ' ' V 

■ , n ■ ....  ..  /,;  •"  i:;, 

i^^The  decimal  point  in  the  coefficients  A**  .Ax'.Ay'  and  A &*!!'  was  moved  for 

i1  O'  .0  . Q 

the  convenience  of  the  numerical  computations.  « 

■ / 1 ' \.  , 't.'  " - . "■  k " 

11  II  » “ I!  ■■  V ■■■'  //  ■ ■ , " ■ „•  \lii  " 


— ‘ - .V  A '".‘A  _ 


the  normal  equations 


jy  ? • v 

, fl  ■ 


A a A®  |i 

Ak  !'  AX„ 

0 

AY0  * 

A 7 

>.  -Q 

♦.2796869  -.00001^3 

+so5H63ii5  -1.2507295 

-.0001203 

. •' 

-,036iiU35 

♦.2785281 

-.0685966 /41  +.000191*1*  +1.21*621*96 

-.0281270 

♦1U2U3535  -.21*70901 

-.32275U9 

+ .000031*5 

ii  ■' 

U *"  /V  ’ • " 

vs  +5.6030266 

0 - 

♦.1518191* 

u )) 

?!  »•//•  •... 


-.1116531*  -.OOii3Ql*l*99Q46 
♦ ,3935762  -.001263766873 
+ .00001391*6912 


Solution  of  the  normal  equations  gives: 


and  from  (?5^ 


Ac  - 4|0107C070  AIL  * >31*566  Ar  « -.01070070 

A®  - / .00356141  AIq  - ♦17.033  As  « -.00356U11 

A k » ♦ .0133231*1  A 2 - +32.129  1 At  - + .0051*1912 

0 . M „ 


V | A K m -t, 01332 3lil 

•h.  . h ■ * . 


2nd  approximations  from  formulas  (66)  and  (73) » 


9 ii 

*1 

;>  ■ *./ 


I.  5. 

V i 


X°  * 581.56  ‘A  \ !r°  - -.01070070  r,(?»  v 1 - r°  » -.99991*275 

°'-.v  '*  V •:  ..  « ?.  * ' 11 

/•  \ II.  I ■ X 


r - 657.03 

o 1 ■ 


Z°  =2632.13 
o 


S°l-  -.00356U11  s °-  Vl  - s°  =+.99999365 


t°  = + .918961*58  t'°-\/i  - t°  - +.39U3U008 


)/■  * .V  '■> 

t.L  f 

■is  11 

- *!  r us./ 


A0 

-+.391*35255 

a’°=  +.91889693 

D°  -'-,01070063 

Vi  j 

B° 

-'-,91895871* 

B*  °=  ’+‘>,391*33758  la 

!l  ‘1'  m ;■ 

E°  - ->053561*11 

C° 

- -.00091*1*61 

c* c*  -,01123696 

'i  ii  " 

•i  •>  * ..  ij 

F°  - 9959361*0 

t % 

f ' /:/!, 

*"  'I  1/  • \v 

. ii  ' -■  £■ 

■r  u 

...  " . // 

■ 7 >r 

[i  • ■<  • ' - //: 

; / ' f)  ■ •*  ■ ■ 

//  ’ 

lj  - ? \ " 

.1  v 

ica* 


l,;' 


l.  i"v\ 
'S 


in 


.u  v 
6 


+557.0521*11 

\ ."■■"!  "!i„ 

+1*81*. 251522 

+129»l*6??f)9 

;\Y  f| 

+i*9 .32511*2 
+63.727711 
-338.675809 
-1*22.057556 
-1*82.527186 


n 

65;.  3067.1  a 

-2,|ai»866l 
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The  normal  equation  are  ; 
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Th©'  weight  coefficients  from!  the,: general  solution  of  the  nomal  equations  are!; 
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The  solution. of  the  normal  equations  gives: 
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mean  error  of  an  observation  of  unit  weight  after  the  adjustment. 


ih\ 


f Confutation  of  re  si  dual. <) 


r 


..  'A 


I « 


Final  Cheek  i! 

■ r '! 

0-(/+v)  o'  -(/Vv'i) 


microns 

,,  " 

microns 

II  :• 

./  8th  decimal  place 

-0,33 

••  1!  " II 

♦1.05  " 

// 

•./<<  ; .7 

/ \\ 

-28 

+6.70 

+5. hit  ' : ' 

, ■ * / -6 

-13 

-13.38 

+1.53 

-|3'  ■'% 

ii 

•?  -9 

•:  ill  : 

. v. 

+12.77 

-0.75  / 

45;, 

0 

1 

-0.1*3 

+19.18 

1)  . *•  V* 

■ ' 

<r\ 

eg 

1 

-2.81 

-9-72 

. 0 -t 

-32 

-1.19  11  X' 

.»  <A\  ;>  '* 

-11.87  "tt 

/J 

■T  +1 

-29 

-7.82 

♦1,60  V 

••  | +10 

V~?5 

+6.1*1* 

-6.26 

^ a • ■ 

, 0 

-20 

[vv]  - 

1178.  , 

’ * ••  '*  •'  .. 

fpLL.6]  - 

■i 

1177.  ji 

?.*ft  'j;,:  // 

* 

L " J 

ii  r • 

- A<  7, 

r.  £ 

■■1 

: ■ ii 


" -i 
X 

■ 

ii 


"il 


/// 
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If  the  prepstablishod  discrimination  constant  ms  * 0.0000001$  the 

*•  iteration  will  be  stopped  at  this  stage.  p- 
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. I ■ -i  ,,  . SOBSAlit  ■ 

,,  A comparison  ipf  the  effort  involved  in  the  dif  ferent  numerical 
solutions  for  the  [orientation  of  a photogrammetric  camera  shews,  that  the 
solution  based  an  the  general  expressions  explicit jjia  the  terms  \pf  the  ob~" 
served  plate  coordinates  (formulas  (h&),  pages  JT-jjh  and  numerical, 
examples 'Nos,  2 and  3)  is  the  host'  economical.  solution  end  suited  for  ! V 
electronic  computing/devices*  This  statement  is  base jl  oh  the  ajqp^piphcd  ■. 
gained  from  a large  number  of  orientations  computed  by  t%  different , ■ 

;; methods*  It  may  be?  seen  directly  by  comparing  the  characteristics'  of  " v.l.yj 
this,  solution  with  /those  of  other  possible  a««roacljes  as  outlined in  ’ ;'V 

tho  present  report  as  well  as  SRI,  Hoeort  Ko*  78b.^  * '4  • f*'  , 

r.  --a  •-* . , ■'  " ■ j 

! ■ (1)  The  formulas  do  not  restrict  the  problem,, to  the?/' protective 

relation  of  ■ftacjj  planes^  but  relate  the  image  plane1! (the  photograph)  to 
the  object [/spade.  ;.l! 

” '.i  ■■  M i ■ •>*''  ' ■ ...  ••  ;;  , l . "I  / # 

(2)  This  solution  has  variables  with  physical  meanings  dterived  , 

from  instrument  readings  (dials  and  levels)  or  independent  camera  call-? 
bration  procedures.  | ' 

..  . • •■■■  ' ' ■ '■  - - . ■ b ;i  y 

(3)  Consequently,  there  are  no  additional  transformations  necessary  • 

| from  plate  constants  to  orientation  elements  and  vice  versa*  '■  [ . ^ 

;•  ••  • ji  ' V\\  » J..  ■ , 

H . \ (U)  Furthermore,  the  solution  may  be  reduced  to  the  actual  number, ^ 

of  unknowns,  thus  reducing  the  numerical  work  to  a minimum.  1 l.?  : 4 

(5)  No  additional  condition  equations  need  be  introduced'.  1. 

1 : (6)  No  additional  weights' are  necessary.  % y 

(?)  The  residuals  v and  v>  are  directly  obtained  from  the  obsorva-  ?V  jj 
tional  equations.  •'  ,;j 

»"  . ■ ' ' „ fy  "...  * f " 

• (8)  The  mean  errors  of  the  orientation  elements  are  directly  ob-  ■ 
■/tain.ed  from  the  general  solution  of  the  normal  .equations,.  ..«• 


I" 


(9)  The  solution  is  suited  for  electronic  computing  devices, 


if 


| ' *■"-  l ' « ■■  -.4  ' . APPENDIX  " '* 

I . ■’?  ■■  . . ';■>/  , . !l 

j After  coding  the  orientation  problem  for  a single  photograiranetric 

camera  based  on  formulas  (61*|  - (75),  the  corresponding  computations  have 
l ", , ' 1 been  successfully  carried  out>on  the  OjfDVAC  ci»tpuifcer'iat  the  Ballistic 
| v;;  Research  Laboratories « A £sx((?Lgl.4i/  iteration  carrying  six, unknowns  for" 

:■  ;i  n » 10  points  takes  35  decor^s.  In  general,  it  will  be! possible  to  ob- 

f J l tain  approximate  values  . fpr  .the  unknowns  which  will  make  the  solution 

j—  " Converge  in  five  iterations  or.  less . . Printing  out  the  result  on  IBM* 

| cards  and  tabulating  about  75  answers  (20  residuals  v and  v1 , 20  A/ 

l and  A/4  values  for  the  final  check,  6 unknowns  and  theif  mean  errors, 

21  weight  coefficients ,„TyyJ  , QIjL.63  and  the  m'san  error;  of  an  ob- 
„ servati.on  of  unit  weight)  takes  another  35  seconds,  Thuis  the  required 

« time  for  the  analytical  solution  .of  the  orientation  problem  for  either  a 

spatial  resection  or  a Camera  calibration  takes  aboiit  3 minutes  or  even 
j .less.  On  the,  following  page  a copy  is  given  of  the  original  print  of 

the  result  of  example  No,  3 as  obtained  from  the  ORDVAC  after  2 iter- 
• | 1 at ions , The  result  shows  complete  agreement  with  the  corresponding  hand 

f--  'J  a."  ;!  computed  values  on  page  83.*  ?.  '.  & ...  * 

■■  «.  : The  IBM-output  is  in  the  'floa ting'd ecimal  point  system.  The  last 

’ „ j two  figures  determine  the  position  of- the  decimal  point  with-  respect,  to 

"T'  ' ,,  the  number- -given  in  the  first  seven  figures . ,j  The  number  of  l&e -last  two 

' ' digits  indicates  how  many  places  the  decimal , (point  must  be  moved.  No 
" „ i . sign  in  front  of  these  last -two  figures  indicates  that  the/ decimal  point 

i - must,  be  moved  to  the  right,  a negative  sign  indicates  a corresponding 

a displacement  to  tjn&  e.«g*  f » • >' 
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